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Abstract 



\ Within of the framework of semiclassical approximation a general formahsm for deriving an ef- 
fective current generating bremsstrahlung of arbitrary number of soft gluons (longitudinal or 



transverse ones) in scattering of higher-energy parton off thermal parton in hot quark-gluon 



' plasma with subsequent extension to two and more scatterers, is obtained. For the case of 
' static color centers an expression for energy loss induced by usual bremsstrahlung of lowest- 
I ■ order with allowance for an effective temperature-induced gluon mass and finite mass of the 
projectile (heavy quark), is derived. The detailed analysis of contribution to radiation energy 
1^ ' loss associated with existence of effective three-gluon vertex induced by hot QCD medium, is 
performed. It is shown that in general, the bremsstrahlung associated with this vertex have no 
^ ' sharp direction (as in the case of usual bremsstrahlung) and therefore here, we can expect an ab- 
d ' sence of suppression effect due to multiple scattering. For the case of two color static scattering 
centers it was shown that the problem of calculation of bremsstrahlung induced by four-gluon 
hard thermal loop (HTL) vertex correction can be reduced to the problem of the calculation of 
bremsstrahlung induced by three-gluon HTL correction. It was shown that for limiting value 
of soft gluon occupation number A^k ~ all higher processes of bremsstrahlung of arbitrary 
number of soft gluons become of the same order in coupling, and the problem of resummation of 
all relevant contributions to radiation energy loss of fast parton, arises. An explicit expression 
for matrix element of two soft gluon bremsstrahlung in small angles approximation is obtained. 
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1 Introduction 



In the third part of our work we complete an analysis of dynamics of boson excitations 
in hot QCD-medium at the soft momentum scale, started in jTJ |2] (to be referred to 
as "Paper I" and "Paper 11" through this text) in the framework of hard thermal loop 
effective theory. Here, we focus our research on the study of the bremsstrahlung of soft 
gluons (transverse and longitudinal ones) by high-energy color particle (parton) induced 
by collisions in the quark-gluon plasma. This color-charged parton can be external one 
with respect to medium or thermal one, and will be denoted by Greek letter a in the 
subsequent discussion. For the sake of simplification we consider QGP confined in an 
unbounded volume. 

Here, we present a formalism permitting to calculate by systematic way the proba- 
bilities of the bremsstrahlung in the semiclassical approximation with taking account of 
thermal effects of medium. These effects are appeared not only in 'dressing' gluon prop- 
agator, but also in 'dressing' of vertices of interaction by thermal corrections. At present 
Gyulassy-Wang model jHlll] used for description of radiation processes in QCD medium is 
generally accepted one. It presents the QGP as a system of Debye screened color centers. 
The Debye screening mass fiD plays the role of a natural infrared cut-off that enables 
us to avoid infrared divergency, which arise when we integrate over transfer momentum 
q. Thus in this case in principal the whole effect of medium is reduced to appearance of 
the only parameter fiD (if forget at the moment about influence of multiple scattering on 
bremsstrahlung). But generally speaking, influence of medium is not restricted by this. 
Fast parton can scatters not only off static color center, but also scatters off thermal par- 
tons which form a screening 'cloud' of static center with subsequent bremsstrahlung of soft 
gluon. In usual plasma such type of the bremsstrahlung was first considered by Akopian 
and Tsytovich jHl Ej and named transition bremsstrahlung. The transition bremsstrahlung 
is a purely collective effect. In hot QCD plasma such type of bremsstrahlung is produced 
in the leading order in coupling (in the framework of semiclassical picture) by precession 
of color classical vectors of thermal partons forming a screening cloud and effective de- 
scribed by the entering of three-gluon hard thermal loop vertex correction ST^g j7j . Thus 
the picture of 'dressed' color particles manifests itself not only in collisions and scattering, 
but also in emission processes. By virtue of the fact that thermal particles of polarization 
radiate, we can expect that angle distribution of such type of the bremsstrahlung have no 
sharp direction as it takes place in usual Compton bremsstrahlung, and therefore it will 
not to be suppressed by multiple scattering - Landau- Pomeranchuk-Migdal (LPM) effect 
|H] . One of the aims of our work is defining frequency range u of emission quanta, where 
HTL vertex corrections can give essential contribution to radiation energy loss. This re- 
gion is restricted by frequencies u of order gT, where T is a temperature of system. In 
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the region uj ^ gT a contribution to energy loss induced by the effective vertex ST^g, is 
suppressed by couphng g. 

Although our consideration is focused on the study of static limit of scatterers, nev- 
ertheless obtained initial expressions for the probabilities of bremsstrahlung enable us in 
principal to extend the results of the work to more general case of moving scatterers. 
At first consideration of such (more physical) situation was made by Arnold, Moore and 
Yaffe [H] starting from somewhat different statements. A motivation to this is the fact that 
typical scatterers are themselves moving at nearly the speed of light, and they produce 
dynamically screened color electric and magnetic fields formed fluctuating background 
field in which bremsstrahlung and the LPM effect take place. Note that a problem of 
importance of accounting of dynamical screening in the magnetic part of the one-gluon 
exchange in the study of radiative energy loss of fast parton propagating through QGP, 
has raised in the work ^Hj- The formalism developed by Arnold, Moore and Yaffe [Oj 
on the background of thermal field theory has valuable advantage in comparison with 
different approaches (see below), since it enables us to correctly treat the bremsstrahlung 
and the LPM effect at all energies of emitted gluon u > gT (as distinct from standard 
requirement u ^ T). In particular it enables us to study an phenomenon of gluon emis- 
sion in a hot QCD plasma induced by collisions of thermal partons between themselves. 
The estimation of energy loss of light partons in the context of approach (01 has made by 
Jeon and Moore p.Tj . 

Besides our approach permits exactly to take into account so-called double Born am- 
plitudes [I2l[l3j. In the works of Zakharov [12] on the basis of path integral approach to 
the induced radiation, it was shown that interaction of fast quark and emitted gluon with 
each color static center will be described by not only one- but double-gluon exchanges. 
The interference between the double-gluon exchanged diagrams and the diagram without 
gluon exchange is important to ensure unitarity. In our approach the contact double 
Born terms arise from 'non-diagonal' elements in correlator of the product of two effec- 
tive currents (Section 5). These terms enable us to cancel exactly divergencies arising in 
main 'diagonal' contributions to energy loss associated with emission of quanta lying off 
mass-shell. 

The main application of approach to bremsstrahlung of soft gluons developed in our 
work, is a problem of energy loss. The research of the energy losses of energetic partons in 
QGP is of great interest with respect to jet quenching phenomenon including in suppres- 
sion of leading hadrons from fragmentation of hard partons due to their strong interaction 
with the hot and dense QCD medium. This phenomenon was long predicted by Gyulassy 
et al. jH! in the framework of QCD-based model calculations and at present it is observed 
at Relativistic Heavy Ion Collider at Brookhaven National Laboratory [15-24]. 
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The first calculations by Gyulassy and Wang |3jj, and Wang, Gyulassy and Pliimer 
jl] of QCD radiative parton energy loss incorporating the LPM destructive interference 
effect indicated that medium induced by non-Abelian bremsstrahlung dominates over the 
elastic energy loss [23 1211 12Z] • This important conclusion of these two fundamental works 
stimulates further more extensive research of mechanisms of radiative parton energy loss 
in hot nuclear matter^. In particular the following important step was made in works of 
Baier et al. [3^, where it was pointed to necessity of taking into account of the effect of 
radiated gluon rescattering induced by multiple scattering in QGP. At present there are 
few developed techniques for the computation of the gluon radiative spectrum and the en- 
ergy loss included a light-cone path integral approach original proposed by Zakharov J2] 
and effective 2D Schrodinger (or diffusion) equation ^33,, opacity expansion framework 
pij (opacity expansion in terms of path integral formulation j35.) supplemented by exact 
Reaction Operator Formalism [SHI which is more suitable for multiple parton scattering in 
a finite system and a twist expansion series [SZ| • The detailed reviews of these techniques 
one can find in EHl HO] • Approximate ways of estimating the observable in- medium 
quenching of inclusive particle spectra include a suppression of the hard partonic cross 
section jH] or an effective attenuation of the fragmentation function 021 03] • One of the 
further developments of the theory of radiative energy loss is connected with consideration 
of more subtle effects caused by the change in the cut-off scale for parton splitting and 
running coupling constant from the vacuum to QGP [l^, the influence of flnite kinematic 
boundaries on the induced gluon radiation pUl I45j . the influence of the color dielectric 
modification of the gluon dispersion relation in QGP (Ter-Mikayelian QCD effect) [IH], 
allowance for mass finiteness of high energy parton ('dead' cone effect) |42i. Let us con- 
sider in greater detail the reviews of the results of works concerned with research of the 
two last above-mentioned effects as more close to subject of present paper. 

The Ter-Mikayelian effect (or as sometimes it calls density effect in the bremsstrahlung) 
for the condensed media |1HI HH] or QED plasma jSl EOj is associated with suppression of 
spectral density of radiation intensity of the bremsstrahlung in frequency range uj -C u;pi7, 
where a;pi is a plasma frequency and 7 = i?/M is a Lorentz factor for particle with energy 
E and mass M. Thus an existence of this effect means not only the regard of interaction of 
emitted quantum with medium (in particular, in gain of medium-induced effective mass by 
it^), but allowance for the finite mass of high energy particle radiated this quantum. For 
hot QCD medium this effect was first considered by Kampfer and Pavlenko [IB] by means 
of the introducing of an effective gluon mass, depending on the temperature of system. 

^Here, we have not discussed the problems associated with calculations of the energy loss of high energy 
partons moving through the cold nuclear matter, which possess proper specific features j28ll2t)ll!MlHlli;i2| . 

^Note that in the work of Galitsky and Yakimets j2] (see, also 122]) it was proved the important fact 
that besides medium-induced mass, the effect of absorption by medium of bremsstrahlung of ultrarela- 
tivistic particles also can result in strong inhibition of bremsstrahlung in a certain frequency range. 
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It was shown that QCD Ter-MikayeUan effect leads to the suppression of the induced 
radiation in a different phase space region, where LPM effect is appeared. Furthermore 
Djordjevic and Gyulassy [SH! have considered a problem of influence dielectric properties 
of hot nuclear medium on the zeroth order in opacity associated radiation and shown that 
the Ter-Mikayelian effect reduces the associated energy loss. The problem of in-medium 
modification of the parton cascade without gluon exchanges between the fast parton and 
thermal partons has also studied by Zakharov [44]. 

Now a study of radiative energy loss of heavy partons (c or h quarks) propagating 
through hot QCD medium is of great independent interest. At first Shuryak [54^ have 
pointed to importance of taking into account of the effect of the energy losses on heavy 
quarks in A + ^ collisions in the context of their influence on high-mass dileptons from 
heavy quark decays. These dileptons can stand out as one of penetrating probes of 
the most dense stages of nuclear collisions IH7j . Furthermore in the work of 

Mustafa, Pal and Srivastava [SB] semiquantity estimation of the radiative energy loss 
of heavy quarks was given, and it was shown that this energy loss is larger than the 
coUisional energy loss [SHI lEOl IMl for all energies of practical importance (see, however 
|68j). Finite mass effect was taken into account in jHE] through kinematic bound on the 
maximum transverse momentum of the emitted gluon without modification of multiplicity 
distribution of the radiated gluon. 

The bremsstrahlung of heavy charged particles possesses proper specific features as- 
sociated not only with different kinematic boundaries in comparison with light charged 
particles. In the case of QED interaction at first a problem of qualitative difference of the 
bremsstrahlung of heavy ultrarelativistic particles (type of muon or proton) in dense an 
amorphous medium was considered by Gurevich It was shown that taking into ac- 
count of mass of projectile is not reduced only to replacement of electron mass by mass of 
proton (for example) in formulas of the bremsstrahlung. The main principle conclusion of 
the work jHlI is consist in the fact that from four main ranges of the bremsstrahlung with 
different frequency dependence of spectral density of radiation intensity ((I) the range of 
absence of medium influence (Bethe-Heitler); (II) the range of medium polarization (Ter- 
Mikayelian); (III) the range of absorption of photons [SI1E2]; (IV) the range of multiple 
scattering (Landau-Pomeranchuk)) for heavy particles multiple scattering range vanishes, 
and three first range of bremsstrahlung remain. At qualitative level this occurrence can 
be explained by that in the case of a heavy particle the formation time of photon radiation 
is reduced to a light one, that in turn results in a significant reduction of the LPM effect. 
This conclusion is also hold for particles with non-Abelian type of interaction. 

In the case of QCD interaction a problem of qualitative difference of radiative heavy 
quark energy loss from that of the light quarks, was considered by Dokshitzer and 
Kharzeev [Hj. Here, it concerned with suppression of small- angle radiation due to so- 
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called 'dead' cone phenomenon ^j. The regard of the finite mass results in appearance 
in the quenching factor |¥71 inE] supplemented positive term specific for heavy quarks, 
that leads to the reduction of the suppression of the heavy hadron p± distribution in 
comparison with one for the light hadrons. The work of Djordjevic and Gyulassy jHH 
was the next important step in the study of the problem of heavy quark energy loss. In 
this paper it was be proposed an extension of opacity expansion for massless partons to 
heavy quarks with mass M. It was shown that a general result (taking into account the 
Ter-Mikayelian plasmon effect for gluons) is obtained by a single universal energy shift of 
all frequencies by Aa; = (m^ + x'^M^)/ {2xE) in the Gyulassy-Levai-Vitev series [SHEE], 
where rrig is an effective plasmon mass and x is a light-cone momentum fraction. The 
importance of developed approach in the work [HZ] is determined by the possibility of its 
subsequent application to development of the theory of heavy quark tomography of hot 
nuclear matter^. 

In the work of Thomas, Kampfer and Soff [ZO] a numerical calculation of exact radi- 
ation amplitude for heavy particle as function of radiation angle 6 for the case of single 
scattering, has performed. It was shown that although there is a suppression effect due 
to the heavy quark mass, its simple consideration by entering the dead cone factor is not 
correct in all kinematical situations. Besides a test of the validity of the potential model 
beyond the assumptions of massless scattering particles by direct calculation of the radia- 
tion amplitude in a quark-quark collision, when a target quark is at rest, was made. The 
results of the calculation shown a departure from those based on Gyulassy- Wang model 
for large angles and when the projectile is heavier then the target, i.e. when it cannot 
be neglected by transfer of the energy to the target, and here, such the necessity of im- 
provement of potential model arises. In the works of Armesto, Salgado and Wiedemann 
|7Tj and W.C. Xiang et al. \T2, similar problems were considered within the path integral 
formalism. 

The paper is organized as follows. In Section 2 a general formalism of deriving an 
effective current generating the bremsstrahlung of arbitrary number of soft gluons for 
scattering of two color partons among themselves in hot QCD medium, is introduced. 
Section 3 is concerned with deriving a formulae of radiation intensity of the bremsstrahlung 
of soft gluon generated by the lowest-order process of the induced gluon radiation. In 
Section 4 the expression for radiation intensity derived in previous Section is analyzed 
in a simple case of the potential model, and under condition when HTL-correction to 
bare three-gluon vertex can be neglected. Section 5 presents a detailed consideration of 
the energy loss connected with bremsstrahlung induced by effective three-gluon vertex 

■^The problem on modified fragmentation functions and radiative energy loss of the heavy quarks 
propagating in a cold nuclear matter was discussed by Zhang, Wang and Wang |68l in the framework 
of the generalized factorization of multiple scattering processes. 
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ST^g. In Section 6 a problem associated with an existence of off-diagonal contributions 
to radiation energy loss is considered, and their connection with so-called double Born 
scattering is discussed. In Sections 7, 8 and 9 an extension of the method of the effective 
currents to high-order processes induced by soft gluon radiation: bremsstrahlung of two 
gluons and soft gluon radiation in the case of two-scattering thermal partons, is presented. 
In Section 7 derivation of an effective current generating bremsstrahlung of two soft gluons 
is presented, and corresponding general expression of intensity radiation in this case is 
written out. In Section 8 an approximate expression in small angles approximation and 
in static limit is obtained from general expression for matrix element of two soft gluon 
bremsstrahlung. Section 9 is concerned with derivation of an expression for gluon radiation 
intensity for the case of scattering of a high-energy incident parton off two hard thermal 
partons moving with constant velocities. The case of bremsstrahlung induced by four- 
gluon HTL vertex correction is analyzed in detail. In Section 10 results of this work is 
briefly summarized, and some ideas of their extension to the case of dense strongly coupled 
quark-gluon plasma actually observed at present at RHIC experiments, are proposed. 



2 Effective currents for induced soft gluon radiation 

One expects the word lines of the hard color particles to obey classical trajectories in the 
manner Wong [73] since their coupling to the soft QCD-plasma modes is weak at high 
temperature. Considering this circumstance, we add color classical currents of hard point 
charges a and 

= KOT)^^'nx-xo„-v,(t-to)), jJ^(X) =^?t;^g;(t)5(3)(x-xo;3-v^(t-to)) 

(2.1) 

to the right-hand side of the field equation (1.3.1). Here, f^^^ = (1, v^^j); Qa,i3 are color 
classical charges of the particles a and f3; xqq^/j are their initial states at time t = to- The 
color charges satisfy the Wong equation (II. 3. 12). Its solution, for example, for particle a 
takes the form 

QUt) = U^%to)Ql, QL = Qam=to, 

where ^ 

U{t,to) = Texp{-zg [ {v^ ■ A'^(r, v„r))T'^rfr}, {T'^f' = z/'^^^ (2.2) 

Jto 

is an evolution operator taking into account the color precession along the parton trajec- 
tory. Rewriting the field equation in momentum space we lead to the nonlinear integral 

^In what follows the high-energy particle a will be considered as external one with respect to medium. 
It either injects into medium or produces inside whereas particle (3 (or particles /3i, f32, ■ ■ ■) is hard thermal 
particle (or particles) with energy of temperature order T. 
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equation for gauge potential instead of Eq. (11.3.2)^ 

*v-'>^'^{k)Ai{k) = -JtrAAm - - fQ^,[A]{k) ^ -j^^^^^^^iAm. (2.3) 

Here, 

00 

^rL[^](fc) = E'^^^^""(A...,A), (2.4) 

s=2 

A'^^A, . . . , A) = i / *T^::Z{k, -k,, . . . , -k,)A^^^^^{k,)A''''Hh) . ■ ■ a^^^^^s) 



s\ 



s=2 

flfj,l...^ls 

s s 

X 6{k -Y,ki)I[dki 

1=1 i=l 

is nonlinear induced color current, where the coefficient functions *rf"/ ""f are usual HTL 
amplitudes, and 

^ """S (2^-/ {va-{h + ... + k.M'^a ■{k2 + ... + k,)) . . . K ■ A;,) ^^-^^ 

s 

x{v^ ■ A^^iki)) ...{v^- A'^^iksMva ■ {k - ^fc,))e'^'"S'''^"" n^A;, (T"^ . . .T^^^r'Q',^ 

i=l i=l 

is Fourier transform of color current of hard particle a ()2.ip . ()2.2|1 (a similar expression 
holds for particle j3). In derivation of two last lines in Eq. ()2.5|) we drop all terms containing 
initial time to- *'D^^^{k) is a medium modified (retarded) gluon propagator which in the 
considered problem for convenience is chosen in Coulomb gauge *T>'^'^{k) = *'D(^{k). In 
the rest of frame of medium we have 

(^) = (^) *^\k)- *V''\k) = 0, (2.6) 

*v'S{k) = {5'^ - uy) *A\k), je = ky\k\, 

where *A''*(A;) = l/{k'^ — n''*(A;)) are scalar longitudinal and transverse propagators. 

As in Paper I (Section 5) and Paper II (Section 3) the first step is consideration of 
a solution of the nonlinear integral equation ()2.3p by the approximation scheme method. 
Discarding the linear and nonlinear terms in v4^(fc) on the right-hand side of Eq. ()2.3|1 . we 
obtain in the first approximation 



*V-'^-^{k)A:ik) = -M^'ik) - 4t^(A:), 



^In the third part of our work we replace a notation of 4-momentum p for Fourier-transformation by 
k. Thus our work in these notations will be close to most of works considered QCD radiation energy loss. 
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where Jq^a^pik) = g / {2t[Y QQ^v'^p5{va^p ■ k)e'^^'^°'^''^ are initial currents of high-energy 
parton a and test thermal parton /?, respectively. The general solution of the last equation 
is 

A^ik) = Af-{k) - *Vc,Ak)J^Sl^'{k) - *Pc..(A:)4t'(^)' (2-7) 

where A^^^'^ik) is a solution of homogeneous equation (a free field), and two last terms on 
the right-hand side represent a gauge field induced by hard partons in medium. 

Furthermore, we keep the term quadratic in field in nonlinear induced current ()2.4|1 
and term linear in field in 'precessional' current ()2.5p . Substituting derived solution ()2.7p 
into the right-hand side of Eq. ()2.3|) . we obtain next correction to the field ()2.7p 

- *Vc,Ak)[iT''TQL^, J K ■ ■ {k - k,))dk, 

-*Vc,Ak)[{T''r'Qlo.j^, I K^-DcMkVt^'^'ikrmv^ ■ (A: - k,))dk. 



2 ,. 



Here, the first term on the right-hand side is associated with a pure soft gluon self- 
interaction. It was analyzed in Paper I. The second group of terms can be rewritten in 
the form 

where effective current is 

jSa'i^^'^m = (T^'r'QLj^Ys Jk^'-m k, -kMTKk,)6{v^ ■ {k - ^0)^^!, (2.9) 

and in turn 

K^^\^r^\ k, -k,) = /"""f . + *T^^^''{k, -k,, -k + ki) "Vc^k - kiXa . (2.10) 

{Va ■ ki) 

The current ()2.9|1 induces so-called nonlinear Landau damping process, studied in detail 
in our previous Paper II. 
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The third group of terms on the right-hand side of Eq. ()2.8|) is equal to zero by virtue 
of identity f°'^'^Qoa isQoa p ~ 0- Now we concentrate on new nontrivial terms, forming the 
last fourth group. After algebraic transformations they can be presented as 

- *1^C^Lu{k)K^{-Va, V^; Xoa, Xo/3, Qoa, Qof^l k). 

Here, for the case of linear dependence on Qoa and Qo/3 we have 

Kli^a, Xoa, Xo/3, Qoa, Qo/3| k) = K^^^^ a, ^fS] XQa, ^0^1 k)QoaQof3^ (2-11) 

K;''^(v„, v^;xoa,Xo;3| k) = - (T")^^fs:^.(v„, v^;xo„,Xo/3| fc), 

where 

-f^/i(Va, V/3;Xo„,Xo/3| k) = /i {Vau*1^C^' {q)Vf3u') JJ^ T {vpu*1^c''' i^^ " q)Vau') 

J iVa ■ q Up ■ \K — q) 

(2.12) 



New effective current ^^'^(vq,, v^; . . . | /c) generates the simplest process of bremsstrahlung 
of soft gluons. On Fig.^ diagrammatic interpretation of three^ terms on the right-hand 



+ f k - q 




Figure 1 : The simplest process of soft-gluon bremsstrahlung generating by color effective current 
()2.1ip . The blob stands for HTL resummation, and the double lines denote hard particles. 



side of Eq. fl2.12|) is presented. 

Thus we have shown that the effective current generating the lowest-order process of 
the induced gluon radiation from fast partons appears in the solution of the basic field 

^In the semiclassical approximation the first and the second terms on the right-hand side of Eq. H2.12|l 
also contain processes, where the soft gluon is emitted prior to the one-gluon exchange graphs which we 
drop on Fig.^ 
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equation ()2.3|) that defines interacting soft-gluon field in the form of an expansion in 
an initial value of color charges Qoa and Qo/3- Based on this result and on the results of 
Paper I and Paper II also we can write out now more general structure of the effective 
current in the form of a functional expansion in a free field A^^^ and color charges Qoa, 
QojS generating the bremsstrahlung of arbitrary number of soft gluons for scattering of 
two hard color partons among themselves 

jQ^^A^^k) = + Jj^llik) + K;{y., v^; Xo„, Xo;,, Qo., Qop\ k) (2.13) 

oo „ n 

+E tai:::^ ^/^^ ^oa, Qoa, Qop\ K -ki, -k)a^''>'^''^^ (h) . . . a^''^'^-^- {K)l[dks. 

n=r s=l 

The functions -ft'^)^^ "^" , n = 0, 1, . . . itself represent in general infinite series in expansion 
over color charges Qoa and Qo/3- Thus, for example, for K^i^ we have 

-^^^(Va, V/3;Xoc,,Xo/3,(5oa,(3o/3| = ^^''(Va, V;3; Xqo, X0/3 | /c)Qo„Qo^, (2.14) 

+ ^«'^'^Hv„, v^;xo„,Xo^| A;)g^,gS;3gS^^+irf ^"^^(v,, v^;xo„,Xo^| A;)QM?.Qo/3} + - • • • 

Here, the first coefficient function i^^'"^(vQ,, V/j; ... | A;) is exphcitly defined by Eqs. ()2.1H) . 
()2.12|1 . and exact form of next-to-leading coefficient functions ii'^*^^'^^(va, v/?; . . .\k) and 
K'l^^^'^'^iyai V/?; • • • I ^), and their physical meaning will be given in Section 6. 



3 Radiation intensity of soft gluon bremsstrahlung 

In the expression for an effective current ()2.13j] without loss of generality one can set xqq = 
0, and choose vector xq/? in the form xo/3 = (b,zo/3); where two-dimensional vector b is 
orthogonal to relative velocity v„— v^. Besides in the subsequent discussion a longitudinal 
component zq^ also does not play any role and thus it can be set equal to zero. The energy 
of radiation field E§*(fc) = -icj *P^^(fc) jg,'^°[A(o)](fc)+ifc* *Pg^(/t) jg,'^"^ generated 
by the effective current, is defined by an expression 

W{h) = -{27^YjdMujdQ,Mp^l^{jt^*\Khyv^J{k)J^^^^^ (3.1) 

Here, 

dA 

dQoa,^ = n dQoa,P 5(QoaQS/3 - Ct'''^), dA = Nj - 1, 
a=l 

with the second order Casimir for a {(3) color particles; (■) is an expectation value 

over the equilibrium ensemble. Eq. p.lj) represents a minimal extension of an correspond- 
ing expression in Abelian theory [50j taking into account a color degree of freedom of 
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particles. The radiation energy ()3.H) is a function of two-dimensional vector b being 
impact parameter in this case. 

Let us define an expression for radiation intensity of soft gluon bremsstrahlung in 
passage of high-energy parton a through medium consisting of hard thermal partons of 
/3 sort, using basic formula ()3.1|) . Here we follow reasoning of Ginsburg and Tsytovich 
|50j . For simplicity we assume that particles /? are non-relativistic, their distribution over 
momentum is given by distribution function /p^. It is defined by formula 

=Jfp0 



(2vr) 



3 ■ 



In this case radiation intensity is read 

We note once again that this expression defines a change of energy field in system for 
per time due to bremsstrahlung process. It equals (with opposite sign) change of kinetic 
energy of a particle and (3 particles in medium on which nonelastic scattering arises. In 
this case when both particle a and particles f5 radiate it is impossible in principle to 
separate a contribution of energy loss of particle a from a contribution of energy loss 
of particles /3 in medium. In limiting case "frozen" f3 particles only, when it can be 
neglected by its bremsstrahlung, Eq. ()3.2|) coincides with an expression for energy loss 
X|v^=o = —dE/dt, where E is energy of a fast parton a. 

One can define a probability of bremsstrahlung. Let us refer the probability to unit in- 
tervals of momentum of emitted quantum (ik/(27r)^ and transferred momentum dq/ (27r)^. 
The probability of bremsstrahlung w;p^,p^(k, q) is defined by a simple formula 

/ dW{h) \ f /p/P/3 dkdq 



Now we obtain an expression for radiation intensity generated by the lowest-order 
process of the induced gluon radiation (Fig.[l}. We substitute an effective current ()2.1H1 
into the right-hand side of Eq. ()3.H1 . Using an expression for propagator in Coulomb 
gauge ()2.6j) and averaging rules over initial values of color charges 

r C^"'' r 

I dQoa QoaQoa = ~3 ' / ^Qop Q013Q0/3 = "1 — ) 
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from Eq. ()3.1|) we find 



(3.4) 

- WV^ ( v) (^) E/'^k,^- i-CA'(*))l A-V„ v.; b| .)l 

Here, = g'^/'iir. The first term on the right-hand side is connected with bremsstrahlung 
of transverse gluon, and a second one is associated with bremsstrahlung of longitudinal 
gluon. For transverse mode we introduce polarization vectors e*(k, (^), ( = 1, 2 possessing 
properties 

k -6(^,0 = 0, e*(k,C)-e(k,C')=5cC'- 

Three-dimensional transverse projector in ()2.6p is associated with polarization vectors by 
relation 

P'^ik) = {6'^ - Uy) = Y: e«(k, C) e^(k, C). (3.5) 

C=l,2 

Let us consider in more detail an expressions |e*i^*| ^ and \K^ \ ^ on the right-hand side 
of Eq. ()3.4p . Due to Eq. ()2.12|) the first expression can be presented in the form 



I e'^\'{k, C)/C*(v„, v^jI k, g)5(^-v^-q-v„- (k-q))dq 



\e'{kX)K\v^,vp-Mk)\' " ■ 



go=V/3-q 



= f e^(q-q')-b [(^.(k, C)/CXv„, k, g))(e^(k, C)/C^(v„, v^| k, q))*] (3.6) 
X (5(u; - V/3 ■ q - ■ (k - q)) (5((vq - V/j) ■ (q - q))dqdq', 

where 

/C^(v,, v^l fc,g) = ^ (t;«/Pr'(g)^;/3.') ^ r {vp,*V'S''{k - q)v^,,) 

Va ■ q Vj3 ■ [K — q) 

+ *r'^-^(fc, -q, -k + q) *Pc..'(g)< *Vcxx'{k - q)v^. (3.7) 

In the last line of Eq. ()3.6|) delta-function 5((va — v^j) ■ (q — q')) = 6{q\\ — q'^y)/\^a — 
V/3I enables the integration over longitudinal component of vector q' = {(i'^,q\\) to be 
performed. For deriving radiation intensity we multiply ()3.4j) by /p^(ip/3/(27r)^|va — v/jl^b 
and integrate over p/3 and b. We perform integration over impact parameter b with the 
help of expression 

J e^(q-q') b^b = (27r)25(2)((q _ q')^). (3.8) 
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The Eq. ()3.8|) enables us to performe complete integration over dq' in ()3.6p and thus to de- 
fine initial for subsequent analysis expression of radiation intensity for the bremsstrahlung 
process dipected on Fig.[T] 



/ dW{h) \ 
\ dt 



dA 



9r d\pi3\ 
P/3/IP/3I 27^2 



(3.9) 



/^E Jdkdu;ulmCA\k)) |rfq |e*(k, 0^ (v„, A;, g)| ^ 



C=l,2 

X - ■ q - v„ ■ (k - q)) 



If there are several sorts of thermal partons, then on the right-hand side of Eq. ()3.9|) it 
should be summarize over (3. In notation ()3.9p we assume isotropy of a medium, i.e. 

Now we consider a special case of Eq. ()3.9j) . Let us define the radiation intensity 
caused by bremsstrahlung of real quantum of oscillations, i.e. oscillations lying on mass- 
shell. For this purpose, for a weak-absorption medium, when Im(*A~^*''(A;)) — * 0, we can 
approximate imaginary part of scalar propagators in the following way 

ZM(k) 



Im(*A*''(A;)) ~ -7rsign(u;) 



t,i 



(3.10) 



where ;(k) are the residues of appropriate scalar propagators *A ' (k) at the poles 



and cu^'' = co'*''(k) are the dispersion relations for transverse and longitudinal modes. The 
term with 6{u! + cj^'') on the right-hand side of Eq. ()3.1()j) takes into consideration not 
emission, but absorption of oscillation. Substituting ()3.10p into ()3.9p and omitting the 
last contribution after integration over du, we find instead of ()3.9p 



/ dW{h) \ 
\ dt ^ 



X 



dQy,^ 

4:71 



E Jdk 



C=l,2 



/3=(?, g 
/ Z,(k) 
V 2^k 



^2 ^2 
dA 



2 ^ d\pp\ 
P0\ 2vr2 



(3.11) 



on— shell 



X 



6{ujI - ■ q - v„ ■ (k - q)) + (t 



/ k'^ 
V k^ 



<?° = (v/3-q) 



Comparing the last expression with ()3.3p we derive an explicit form for the probabilities 
of soft-gluon bremsstrahlung 



<,p^(k,q) = 4(27r)^a3^Al 



d. 



|e^(k,C)/CXv.,v^|fc,g)|; 

C=l,2 



(3.12) 
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<„/k. q) = 4(2.)«af (g) |/C«(v„. v,| t. ,)| 

The expressions (jH.lip . ()3.12|1 were obtained with the assumption that QGP represents 
a system of non-relativistic thermal partons. In this case only the use of the impact 
parameter and averaging over it, is valid. It enables us to diagonalize the product of 
amplitude e*/C* and complex conjugate amplitude (Eq. ()3.6|) ) in the q, q' variables, and to 
write the probability in the simple form ()3.12|) . In the work [20] a possible way of extension 
of this approach to the case of arbitrary moving thermal particles, was suggested. It 
consists in construction of subsequent Lorentz transformations for particles group with 
momenta in interval (p/3, pp + dp^) at rest, where Eq. ()3.1H) is true. However practical use 
of this idea is very difficult. In work by Akopian and Tsytovich |3] for the case when both 
particle a and particles /5 are ultra-relativistic, the probability of bremsstrahlung ()3.12|1 
for ordinary plasma was derived by direct computations with use of dynamical equation 
without notion of impact parameter. It gives some ground to suppose that expressions 
()3.11|) . ()3.12|) hold for ultra-relativistic high-temperature plasma also. 



4 Approximation of static color center 

Let us analyze an expression for gluon radiation intensity ()3.1H) in the case when target 
particle is modeled by static screened potential^ (v/j = 0), and HTL-correction ST^g to 
bare three-gluon vertex can be neglected. For simplicity we restrict our consideration 
to the case of transverse gluon radiation. At first we consider integral over momentum 
transfer q on the right-hand side of Eq. p.lip presented as 

[dq^dqii e'{k,C)}C\^^a,0\k,q)^ , - • k + w^gy), (4.1) 

where q_L and q\\ are transverse and longitudinal components of momentum transfer with 
respect to velocity Vq,, correspondingly, and Va = |vq,|. Based on an explicit expression 
for /C*(vq,,v^| k,q) (Eq. ()3.7p ) and limiting static expression for propagator 

*^^c(0,q) = --^, (4.2) 

where is Debay screening mass, and performing integration over dq\\, we find instead 
of dHJ 

^ ^ -2 (4.3) 

fqi + /^D + ^(4-va-k)2j 



^Remind that in this case H3.11|l correct to a sign, coincides with the expression for energy loss of 
high-energy parton a. 
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e • Vn 



- ■ k 



90=0, g||=-^(a;^-Va-k) 



Here, for brevity we use a notation e = e(k, (). We note that in our approach in deciding 
on Coulomb gauge, a contribution of thermal parton /3 (the second term in function /C*, 
Eq. ()3.7|) ) to bremsstrahlung of transverse soft gluon is exactly equal to zero. However 
here, there is a possibility of radiation of longitudinal plasmon from target (3. Let us 
consider a difference cj^ — ■ k. We approximate a spectrum of transverse oscillations 
a;^ in the limit of |kx| <^ k\\ by standard expression 

k^ 



- + m2 - kii + -^j^- 

Here, is an effective gluon mass square, depending on temperature of the system, 
= 3u!pi/2, where u^i = g^T'^^Nj + 2Nc)/lS is plasma frequency square. 

Besides, for ultrarelativistic particle a with energy E and finite mass M, we have 

M2 



Based on above-mentioned we can write 



(4.4) 



4-v.-k:^^ , (4.5) 

Furthermore from a condition of transversity we obtain e\\ ^ — (k_L ■ g^_)/uj whence it 
follows 

e ■ v« ^ -Vo, . (4.6) 

UJ 

Taking into account ()4.5j) and ()4.6j) we can write a first term in expression under module 
squared in ()4.H|1 as 

e ■ v„ k I ■ e I , , 

- ~ —2v - - (A 1\ 

oj^-v^-k" " ki + m2 + x2M2 ' ^ > 

and Coulomb factor reads 

1 1 



Here, If = VaTj, where Tf = 2u/{]s.\ + + x^M^) is finite formation time with regard 
to temperature induced gluon mass and finite mass of projectile a. 

Let us consider now a contribution to matrix element associated with bare three-gluon 
vertex. Taking into account structure of propagator in the Coulomb gauge (Eq. ()2.6|l ) and 
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approximation ()4.6p . we can present this contribution as 



-2 



K-(k-q)n 
(k - q)2 



{e-q)*A\k-q) + 2uj{ -v^ 



k_L ■ e_L 



UJ 



e • q 

k-q) 



*A*(A;-g). (4.9) 



Furthermore by virtue of ()4.5|) we have 



kx 1 



e ■ q = e_L ■ q± + e\\q\\ ~ e_L ■ q± + — — j ^ ■ q^ 

under condition Ik^l/cj/j ^ |q^|. Besides we have (k — q)^ = (k — q)5_ + tu^/f^ — uj^/v^, 
when a;^/f^ ^ (k — q)5^. By this mean expression ()4.9p can be replaced by approximate 
one 



x^M^ + (k - q)i 



e± ■ q±) *A\k -q) + 2v^{e^ ■ {v^q - k)^) *A\k - q) 



Here, a first term connected with a contribution from longitudinal virtual oscillation is 
suppressed with respect to a second term by small factor and it can be omitted. 

Let us consider transverse gluon propagator *A*(co',k — q). Its explicit expression in 
hard thermal loop approximation is 



*A*(a;,k-q) 
where 



" -(k-q) -2^^|(k3^-- 



(k-q) 



In 



l + z 



1 - z 



(k - q)2 

-Z7r^(l- 1^1) . 



|k — q| 



By making use approximations 

(k - q)2 ~ k^ + (k - q)i + + x^M^ 
we derive for transverse propagator 

*A*(cu,k-q)~- 



UJ 



|k-q| 



[(k - q)i + x^M^] 
2k2 



-1 



(4.10) 



, ^2 2,.2 1 2 f. [(k-q)i+a;2M2] / 4k2 
k-q)i+.^M^ + -;.^l + ^^i^L^ ^(l^KkZ^i)lT^ 



-ZTT 



The factor before a logarithm on the right-hand side of the last expression is a small 
owing to k^ ^ (k — q)^ + x^M^. On the other hand, the logarithm can be large by the 
same condition. Therefore this term can be discarded under more severe constraint for 
the magnitude of momentum |k|: 

elne<l, e = [(k - q)i + x^M^l/kl 
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Under the last condition with allowance for yu|)/2 = we derive finally expression for 
transverse gluon propagator 

•^''"•■--q'"- (k-q)i+',ng + x^Ap - ("1' 

By this means a contribution to matrix element connected with bare three-gluon vertex 
in this high-frequency limit can be presented as 

eT^^^ik, -g, -k + qyVcAk ' - -2.. (i, ^7+, "2^2 " (4.12) 

In this approximation the residue ^t(k) can be replaced by unit. 

With regard to ()4.7j) . ()4.8|1 and ()4.12j) the expression for energy loss of parton a, 
associated with bremsstrahlung of soft transverse gluon, takes following form within the 
potential model 

^ = 2^f^y^"^^A E E /^k^/^^ (4-13) 



dt J \ n 



e±(k,C) 




The expression fl4.1H|l takes into account the QCD Ter-Mikaelian effect associated with 
existence of an effective gluon mass m^, depending on the temperature T. As was men- 
tioned in Introduction for hot QCD plasma this effect was first considered by Kampfer and 
Pavlenko [IB]. As distinct from this work, in our paper the effective gluon mass appears 
not only in the first term of matrix element, but also in the second one corresponding to 
the radiation from gluon line via the triple vertex. Besides this expression takes into ac- 
count mass finiteness of the projectile (mass effect) in according with expression obtained 
by Djordjevic and Gyulassy [HZl for energy loss of heavy quarks. 

Let us consider more closely the expression ()4.13j) . First of all we define conditions 
wherein a term /^^ in Coulomb factor 1 / (qi+/i|)+^j ^)^ can be neglected. The requirement 
iAd ^ ^/ ^ results in 

k2 1 

-^</iD, -AiD<^<27VD, (4.14) 
luj 4 

where we enter Lorentz-factor 7 = E /M. Performing summation over polarization state 
of radiated transverse gluon, we rewrite the Eq. ()4.13j) as integral over frequency uj 



"dt . 




jdu[r^{uj)+%{u)+%{u)}, 
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where 



r2(cu) 



dk 



dq± 



(qi + /if. + Z/-')^ (k^ 

-2 



kx-(k-q)x 



(qi + /^?) + ^7')^ (ki + m^)((k - q)i + ml) ' 
1 (k-q)i 



(4.15) 



i^l + f^l + lfr ((k-q)i + 



Here, for brevity we denote = + ci;^7 ^. Now we consider a function Ti{uj) 



In 



conditions of ()4.14j) a simple calculation of two-dimensional integrals / dk±J dq^± leads to 



TT 



In 



-I- 777^ 

- max ' g 



1 



(4.16) 



We introduce cut-off k^^^^^^^ on upper limit for integration over dk^. As will be shown 
below the divergent term in a sum of three expressions ()4.15p . is exactly cancelled. Here 
we draw attention to the following circumstance. If we keep the term Ij'^ in Coulomb 
factor, then in this case integrating leads to finite expression by convergence of integral 
over (i|k^|. For conditions ()4.14|1 we will have instead of ()4.16j) 



TT 



In 



TT 

arctan 

,2 



In particular in a frequencies region fi£)/4: <^ u <^ •ynig it follows^ 




In this region a dependence on Lorentz-factor vanishes. 
Finally, in asymptotic region uj ^ 27^/i/) we have 



TT 



1/27Vd^ 



UJ 



Here, we note a strong suppression of spectral density. 

Such allowance for finiteness (inverse) formation length enables us not only to perform 
more accurate integration providing its finiteness, but gives also a possibility to analyze 



^Note that precisely in this frequencies region Ter-Mikaehan effect in usual plasma is manifested. 
It is associated with appearance of suppression factor uj"^ /{uj^ +j^ni'^) ~ {uj/jnig)'^ <C 1 before logarithm 
in spectral radiation density. In the case of QCD plasma such suppression is not observed. 
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a behavior of spectral density outside of framework of restriction ()4.14|) . However here, 
a problem connected with necessity of the keeping of terms of higher order smallness 
over k^^, q^^ etc arises in an approximation scheme, that results in Eq. ()4.13|) . Besides for 
allowance of finiteness of /J"^ for deriving of the remaining expressions of T2{uj) and Tsico), 
calculation complexities arise. In conditions of requirements fl4.14j) before integrating, 
these formulas can be introduced in the form 



-2^ In 



max + "^g 

m2 



TT 



>D - 3m; 



--to 2 \ 

ll In 



ml 



\ml - /i^ 



+ 



2m^ 



In 



fhg + hd 



k^ (ik^ 



(4.17) 



where 



k I max + „iax + TT 



m2 



1^5 -/^D I 



dlC{uj,ml,fil,) 



dw?g 



1 - 



-In 



6 = 2(mJ 



ml 



2 I 2-2 

m„ + u; 7 , 



2-2 2 



D = kl + bki + C, 

JC{uj,ml,fil,) = -^In 

As we see from Eqs. fl4.16|) and ()4.17|) in hmit k^^^^^ — oo cut-off disappears in a sum 
Ti{u) + T2{u!) + %{u) as expected. In region of -C •ymg from exphcit analytic expressions 
there is no any suppression by factor u/'-fmg. 



^In 


rrig + jjD 


CO 


mg 


rhg - 






5 Bremsstrahlung induced by effective vertex correc- 
tion 6r^g 

Let us consider now a contribution to radiation energy loss related to existence of medium 
induced HTL-correction ST^g to bare three gluon vertex. According to Eq. ()3.7p amplitude 
of HTL-induced bremsstrahlung of transverse soft gluon takes a form 

e^(k, C)5r^'(fc, -g, -k + q) *VcuAqH *1^cxx'{k - g)t;^'U=.^,,o=v,.q ■ (5.1) 

In what follows we restrict our consideration to approximation of static color target, i.e. 
we set = 0. As it will be shown below the amplitude 1)5.11) gives a main contribution 
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to energy loss. However it is necessary to note that generally speaking, approximation of 
static color target is in contradiction with conception of hard thermal loops. As known 
HTL amplitude is calculated within high-temperature QCD plasma, when it is assumed 
that all thermal partons (3 constituting plasma, are ultra- relativistic (and massless). Nev- 
ertheless in this Section we consider the medium-induced bremsstrahlung in the framework 
of the potential model, to compare with usual bremsstrahlung within the same model, 
keeping in mind above-mentioned remark. Now we turn to analysis of amplitude ()5.H) . 
In the expression (jS.lj) in the propagator *'Dc\y{k — q) we keep only transverse part, that 
gives also leading contribution in medium-induced bremsstrahlung (see below). Then 
()5.H) in statistic limit reads 



1 



*A*(fc-g)e^(k,C) 



q2 + /il) 



90=0 



(5.2) 



The explicit expression for HTL-correction for go = is 



5r'^'{k,-q,-k + q)\ 



=0 = 1^0^ 



Ait {v ■ k){v ■ {k — q) + ie) 



+0. 



Following Frenkel and Taylor (73] we present integral on the right-hand side as an expan- 
sion in basis 



/ 



V V- 



An {v ■ k){v ■ {k — q) + ie) 



(5.3) 



where n = [q, k] = [k, k — q] , 1 = [n, k] and 1' = [n, k — q] . The coefficient functions in 
this expansion in static limit go = according to Eqs. (3.32), (3.33) in [Tl], have a form 



X = — 



1 

52 



-A M{k, A; - g) + c<;(k • q)L(A;) - cu((k - q) ■ q)L(fc - g) 



Y = 
Y' 

Z -- 



1 

1 
1 

^ 52 



-uj^M{k,k - q) - iul 1 + 



-uj^M{k,k -q) -uj[ 1 



^ k.(k-q) 



k2 

k (k-q) 
(k - q)2 



\L{k) 



k.(k-q) 



\L{k-q) 



k2 

k-(k-q) 



(5.4) 



(k-q) 



-cj^ M{k, k-q)+ ujL{k) + ujL{k - g) - 1 



Here, 



5 = n2 

m = 



A = A;2g2 _ (A; . g)2 = -c^W + 5 



2|k| 



In 







in 


uo — 




2|k 



e I 



UJ 



(5.5) 
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M{k,k 



1 



In 



( 



k-{k-q) + 



k-{k-q)- 



) 



Even in static limit of amplitude ()5.2|) analysis of energy loss induced by medium- 
induced bremsstrahlung in closed analytical form is impossible by virtue of awkwardness 
of the expressions ()5.4|1 . Therefore we make some simplifying assumption for sub- 

sequent analysis. Note that coefficient functions ()5.4|1 contain a factor 1/5^ as a common 
one. Let us assume that main contribution to HTL-induced bremsstrahlung follows from 
kinematic region of momentum variables k and q, where 5 vanishes, i.e. when vectors 
k and q are collinearity. We approximate all functions entering into amplitude ()5.2|1 in 
kinematical region 5 = (5(k, q) = 0. 

In the amplitude ()5.2j) in three-dimensional transverse projector, there is also con- 
tracting HTL-correction with (k — q)''(k — q)'''/(k — q)^. By virtue of Ward identity for 
HTL amplitudes [TH [7S] here, we have 



where 511^'^ is polarization tensor. The vertex 5T'^^^{k, —q, —k + q) on the right-hand 
side contains in coefficient functions the factor 1/6 only in the first power (Eq. (3.30) 
in and such is "less singularity", then vertex 5r*°-^(fc, —q, —k + q). Therefore this 
contribution can be omitted. For the same reason in propagator *T>c\\'{k—q) we neglected 
by the term containing longitudinal scalar propagator *A''(k — q) and in matrix element we 
dropped medium-induced bremsstrahlung of longitudinal oscillations quantum (plasmon) . 
All these contributions are proportional either zero or first power of 1/6. 

We introduce the coordinate system in which axis OZ is aligned with the velocity Vq, 
(Fig.|21). It is convenient to enter a new variable q' = k — q = (q'j_, gy). The longitudinal 
component gy by virtue of delta-function in integrand 1)3.111) (for = 0) equals u/va- 
In future discussion we will write u instead of cj^. In chosen coordinate system we have 
k = (|k|, 6^, q'_i_ = (|q'_i_|, 7r/2, /?), and integration measures are 



Subsequently computation azimuth angles (p and /3 will be enter as a difference (p — (3. 
Therefore we remove integration over j3: J djS = 2tt, replacing if — (3 ^p. The kinematic 
region of variables k and q that is of our interest, is defined by equation 



5r°^(A;, -q,-k + q)ik-q) 



■uj 5r°°(A;, -q, -k + q) + OT°(g) - 6n'%k) 



dk = k'^d\k\ sin 9d9dip, dq'j_ = \q'j_\d\q^j_\dp . 



It is not difficult to see that in this coordinate system one have 





(k-q')' = 0. 



(5.6) 



22 



z 




Figure 2: The coordinate system for analysis of HTL-induced bremsstrahlung. 



— + q_i_ — — coso' + |q^| smt^cosv^ 



J Wn 



0. 



(5.7) 



and using a scalar product defined above, we have 

We consider fj5.7|) as equation with respect to variable cos ip. Eq. ()5.7|) defines two solutions 

cosif± = ?— (coseT Jl + X^), (5.8) 



where we have introduced a new variable x = Va\(i'±\/uj. The condition cos^ V^i < 1 results 
in restriction 



(l^cos^y^l + X^) <0. 



It is evident that the restriction is valid only for two relations between polar angle 6 and 
variable x 



cos 6, 



1 



cos 6'. 



1 



(5.9) 



Substituting these relations into ()5.8j) . we obtain that condition 5 = is true only 
in "end" points: cos^ 99-1- = 1. For definiteness we choose cosip+ = 1, i.e. ip+ = 0, and 
cos(y9_ = —1, i.e. if - = 7T. It follows that sin6'-|- = x/ + X^? ^ind therefore by virtue of 
()5.9p we have 

tan6'± = ±x- 
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Figure 3: The dependence of angles 9± on x- 



The graphs of functions 9± = 9±{x) ci-re given on Fig.El Such the equation ()5.6|) defines two 
non-overlapping cones of medium- induced bremsstrahlung: (1) along a particle velocity 
Vq, (0 < 6-^- < 7r/2) and (2) in opposite direction (7r/2 < 0_ < vr). It is clear that 
such separation exists only within assumption on leading contribution to medium-induced 
bremsstrahlung discussed above. 



Now we define module squared of amplitude ()5.2|) in region ()5.6p . Taking into account 
an expansion ()5.3|) and summing over polarization state of radiated quantum with relation 

E |e(k,C)-a|2 = [k,a]2 

C=l,2 

valid for arbitrary vector a, after simple algebraic transformations we have 
^ •*A\k-q)\' J2 \e\kX)Sr''{k,-q,-k + q)vi\^ 

C=l,2 



where 

J2 |e*5r°^<|' = 5 



C=l,2 



v„-n)2|Xp + (v„-l)2 |k|r+(k-q')Z +(v„-l')' |k|Z + (k-q')r 



+ 2 (v, ■ l)(v, ■ 1') Re({|k|y + (k ■ q')Z}{\k\Z + (k ■ q')r'} j] . (5.10) 

In chosen coordinate system scalar productions of velocity with vectors n, 1 and 1' 
entering into the right-hand side of Eq. ()5.10|) . are represented in the form 

Vq, ■ n = Ci;|k|xsin^^sin(y9, Vq, ■ 1 = cuk^ sin cos [ tan — X cos , (5-11) 
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Vq • 1' = — cjlkl |q'j_| cos6' [x — tan6'cosv9]. 

They vanish for cos^p± = ±1, tan6'-|- = ±x- By using exphcit expressions for coefficient 
functions ()5.4j) one can derive hmiting values of its quadratic combinations which enter into 
the right-hand side of Eq. Rather cumbersome calculations lead to the following 

result 
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where 
M±{uJ, |k|,0 = In^ 
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|kU) - vr^ 



UJ + 


\k\ 


UJ — 


\k\ 



, ^ = lnx + Vl + xM, X = sinh^ 



For convenience of the subsequent reference we write out once more the expression for 
energy loss of particle a connected with medium-induced bremsstrahlung of soft transverse 
gluon 
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X /|q±Mlq±l 



Y: I e^(k, 05r'Kk, -g, -k + g)<| ^ | *A*(4, ql, gf|)|5=^./„^ 

C=l,2 " 



Here, an expression Z)c=i,2 | e*5r*°-'t>^| ^ in limit 5 — 0, is defined by equations ()5.1()j) 
()5.12|) . Coulomb factor in integrand ()5.13j) have limiting value 
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From the equations ()5.10p - ()5.12|) we see that there are four types of integrals over azimuth 
angle 9?, having singular integrands in limit (5 — >■ 0: 



27r 9 2-K 

sm if dip 
53 ' 





2lT 



[tan 6 — X cos p] [x — tan 6 cos p] dp 
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(5.14) 
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[tanO — Xcosp]'^ dp J [x — tan 6^ cos ^9]^ 
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Let us consider for example the first integral. By using an explicit expression for 
function 6 (the left-hand side of Eq. ()5.7|) ) the factor 1/6^ can be presented as 
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P 
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u\k\J (2 cosh 0^ 
+ 3 



(a+ + cos<y5)3 
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+ 3 
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(a+ + 6+ cos<y9)2(a_ + &_ cosp) 
1 1 



(a+ + 6+ cos p) (a_ + 6_ cos pY (a_ + 6_ cos p)^ J 
where a± = cosh^ =1= cos^^, b± = =i=sinh^ sin 6'. The main contribution is connected with 
terms l/(a± + b± cosp)^. Using table integrals [ZH!, we derive 

27r 



sin pdp 



TT 



IT 



[a± + b± cos pY 



(cosh ^ =p cos Oy — sinh^ ^ sin^ ^ 



3/2 



cosh ^ cos 6* =F 1 



3 • 



The integral diverges for values: cos^^± = ± 1/ cosh^. 



Furthermore we consider for concreteness the energy losses associated with bremsstrah- 
lung to hemisphere along velocity Vq, (0 < 6 < '/r/2), i.e. in previous expressions we 
choose upper sign and besides for the sake of simplification we restrict our consideration 
to contribution connected with function X only. By virtue of ()5.10p this contribution 
enters by independent fashion in the squared amplitude, i.e. it doesn't interfere with 
other contributions: Y, Y' and Z. In this case taking into account above-mentioned, we 
have instead of ()5.1()|1 



a;2|k|4 



sinh'^sin'^ |k| coshn MX{uj, |kU) 
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TT 



(2 cosh ^ COS ^)' 



n 3 • 



cosh ^ — 



COS 6* 



By singular factor l/[cosh^ — l/cos6']'^ for rough estimation of energy loss we set 
cosh^ = 1/cos^ in integrand of Eq. ()5.13|) (or in terms of |q'_i_| : |q'_i_| = {oj /va)taxiO). 
The integration measure over |q'j_| will be presented as 

/ 1 q±l d\c{^\ = (-)' f sinh^cie 

J V f rv / COS J 
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We perform the integration of singular factor by the principal-value prescription 

/ l/cos 9—e 



sinh^c?^ 
o = lim 

(cosh ^ - l/cos 
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— 1/ cos( 
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COS^^. 



In the context of this approximation the expression for energy loss ()5.1H|1 associated with 
contribution of function X takes a form 
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and a transverse propagator is 
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Here, the function F(z) was defined in Section 4. 



Let us estimate the expression ()5.15j) in the long-wave limit |k| ^0, c^k ~ ^pi ^ 
small-angle approximation 6 < 7"^, where 7 = E/M. Taking into account 1 — f^cos^ ~ 
(7"^ + 6'^)/2 and integrating over angle 6 in range (0,7""'^), we derive from ()5.15|) 
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^(a)^(/3)- 
^2 ^2 

dA 
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(5.16) 



X (7-2)3 [ln(7-V4) - CJ [21n2 - 5/61 + - In2(7-V4) + 3 ^ 



- (3 + k)k j ' 

where Ce is an Euler's constant and ip{z) = (lnr(2;))', T{z) is a gamma-function. The 
integral over d\k\ infrared diverges in this approximation. It should be cut-off in lower 
limit on ultrasoft scale g^T. It follows 



d\k\ 
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With allowance for the last expression and the fact that for ultrarelativistic plasma den- 
sities rif^ ~ T^, we derive from Eq. ()5.16|) that energy loss connected with bremsstrahlung 
of soft transverse gluons (cu^ ~ gT) in small-angle approximation in value order is 




7 



(5.17) 



e<7- 



Although here, we have the same order in coupling constant as for usual bremsstrahlung, 
however derived expression is strongly suppressed by Lorentz-factor 7. Such in a small- 
angle approximation the energy loss connected with existence of an effective vertex ST^g 
induced by medium, gives negligible small contribution in comparison with usual energy 
loss mechanism for high-energy parton. The most of energy loss, as we expected, is defined 
by radiation of soft gluons for large angles 9 ^ 7^^. 

Now we consider an opposite case of short-wave limit, |k| ^ ujpi, when cj^ ~ |k| + 



m^/2|k|. In approach of small angles and Va 
are obeyed 




cu^ — |k|fQ, cos^^ - |k| 
The Coulomb factor in this case is 



k2 



^1 — 7 the following approximations 
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After integrating over angle 6 in range (0,7 ^) we will have from Eq. ()5.15|1 
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The remaining integrals over (i|k| can be estimated by formula 
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Such in the higher-frequency approximation uj ^ gT we have besides similar suppression 
by Lorentz-factor as in expression ()5.17|1 . an suppression by one power of g. This is a 
consequence of the fact that if one of lines incoming (or outgoing) in HTL vertex ST^g 
is hard, in our case this is external leg of radiation gluon, then this effective vertex will 
be usual perturbative correction to bare three-gluon vertex. This conclusion is certainly 
correct also for energy loss connected with large angle bremsstrahlung. 
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6 Off-diagonal contribution to radiation energy loss. 
Connection with double Born scattering 



The previous Sections were concerned with analysis of radiation intensity of soft-gluon 
bremsstrahlung generated the lowest-order process induced by effective current ()2.1H) . 
It is associated with 'diagonal' contribution ( Jq^^"(/c, b)*r'^'^(A;) Jqj|"(/c, b)) to radiation 
field energy W(h) (Eq. (jH.lj) ). where 

J^^l^k, b) ^ Kf%^r^, v^; b| k)QlQ^,^, 

and function K'jf'^ is given by Eqs. ()2.1H1 . ()2.12|1 . In this Section we briefly analyze a role 
of simplest 'off-diagonal' terms^ 

{j^f^ik, hyv^j{k)jj^l^{k, b)) + {r^f^ik, hrv^\k)j^^l\k, b)) , (6.1) 

where 

4t b) = ^ Ql^v^,5{v^ ■k) + j^, Ql,v,,5{vp ■ k)e^ '> (6.2) 

is initial "bare" color effective current, and jQ^°'(A;,b) is effective current of the next 
higher-order in coupling constant by comparison with Jq^"(/c, b) in the considered problem 
of bremsstrahlung. 

Let us clear up first of all what effective currents of the second order (and the processes 
associated with them) can give nontrivial off-diagonal contribution to energy of field 
radiation p.ip . We consider, for example, a term in the expansion ()2.13p linear in a free 
field 

''[Qoa,Qo/3,A(°)](A;) = y'ir;^i(v„,v^;xo<„xo^,Qoa,Qo/3| A:,-A;i)AW'^^'^Hfci)c^fci- (6-3) 

As will be shown in next Section, this effective current generates bremsstrahlung of two 
gluons in scattering of particle a off thermal parton (3. In the leading order in coupling 
in the expansion of coefficient function i^^^^^ in initial values of color charges we keep 
term linear in Qoa and Q^p only, (see next Section, Eq. ()7.2|1 ). Substituting ()6.3|1 into 
()6.H) and then into (jSH}, taking into account ()7.2j) . after averaging over initial values of 
color charges, we obtain that by virtue of equality / dQoa,i3Qoa,i3 = the off-diagonal 
contribution to the radiation field vanishes. 

Furthermore one can consider an effective current generating bremsstrahlung process 
of soft gluon in scattering of projectile parton a off two thermal partons Pi and l32- This 

^The 'diagonal' and 'ofF-diagonal' terms are defined with respect to those of product of two effective 
currents expansion of H2.13|l type. 
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process will be considered in Section 9. The effective current here, have (in leading order) 
a color structure given in ()9.3p . If one substitutes this effective current instead of Jg^ 
with initial one ()6.2|1 (the last one should be supplemented by initial current of second 
thermal parton) into (jfi.lj) and then into ()3.H1 . at that time after averaging over color 
charges we obtain that this off-diagonal contribution will be also vanishing. 

The only non-trivial 'off-diagonal' contribution to radiation field energy arises from 
the terms of higher-order in color charges Qoa cind Qo/? in the expansion of color effec- 
tive current /('^(vq, v^, ...,Qoa,Qo(3\ k). In this case it is second and third terms on the 
right-hand side of Eq. ()2.14j) . These terms define "classical" soft one- loop corrections to 
bremsstrahlung depicted in Fig.^ We already faced with corrections of this type in re- 
search of the scattering processes of Compton type studied in our previous Paper II. Let 
us substitute the corrections (non-linear in color charges) and initial current ()6.2|) into 
(jfi.lj) and then into (jH.lj) . Performing average over color charges we lead to the expression 
for off-diagonal contribution to energy of radiation field 

(W^(b)):ff_aag. = (^^^j E^yrfkdu;(e(k,C) ■ v.)c.Im(*A*(fc)) 

if"'^^^-'(v<„ v^;b| kyiiiXMv^ ■ k) (6.4) 




d. 



Y: /rfkc/cu(e(k,C)-v^)u;Im(*A*(fc)) 



C=l,2' 



X Re 



e-^'^-^ir'^''"^(v,, v^;b| k)eii}iXMvp ■ k). 



Here, we keep a contribution of transverse part of propagator *'D^'' (k) only and take 
into account expansion ()3.5j) . If the function K"-^^"'^ {y a-, V/3, b| fc)e-'(k, ^) is not singular for 
= 0, then within potential model the second term on the right-hand side of Eq. ()6.4|) 
vanishes. Leaving the calculation details of higher coefficient functions /^^'"=iC2 g^j^^l ]^abib2c^ 
which are similar ones considered in Paper I and Paper II, we give at once an explicit 
form of the first of them 

<^^^Hv.,v^,b|fc) = (|^{^^^^n"''/{~ t;/(r-gi) *'^c'ik-q^)}CA^a,^,\k-quq2) 

+ 7 r^Hw V '(91)^/3.') {VaX*V'c%2M 

[Va ■ [qi + q2))[vcx ■ q2) 

""^^ {vp,*V'S'\qi)vp,') {vpx*Vl}\k - gi - q2)v^x') (6.5) 



{vp ■ {k - q2)){vf^ -{k-qi- ga)) 

+ *r^i.A(fc, -gi, -k + qi) *V''c'{qi)vp^, *Vc^\k - qi)K,x'{^a, k - gi, ga) 

- \ T^,A.(/^, -92, -A^ + gi + g2, -gi) *V-c''{q2)vp,> *V^^'ik -q,- q2)v^x' *Vfiqi)v(,^'^ 
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X e-'^^^+^'>''6{vf, ■ q,)6{vp ■ q2)6{v^ ■ {k - q, - q2))dq,dq2, 

where function /C^(vq,, V/3I k,q) is defined by Eq. ()3.7|1 . Note that integrand is not sym- 
metric with respect to replacement qi # q'2- The diagrammatic interpretation for the first 
term on the right-hand side of Eq. ()6.5p is depicted in Fig. HI 



a ^ 




Figure 4: Some soft one-loop corrections to bremsstrahlung depicted in Fig.0 

We consider further approximation of static color centers, i.e. we set = 0. We 
substitute an expression ()6.4|) with partial coefficient function ()6.5p into formula for radi- 
ation intensity ()3.2j) . Delta-functions in integrands ()6.4|1 and ()6.5j) in static limit results 
in integration measure of the form 

5{Va ■ k)[5{qio)5{q2o)dqiodq2o]5{-Va ■ (Qi + q2))dq,idq2dkduj. (6.6) 

Integration over impact parameter b leads to additional delta-function in integrand 

Jdhe'^^^^-'^-y'' = (27r)2(5(2)((qi + qs)^). 

This expression together with ()6.6j) enables us to integrate over dci2- Performing re- 
placement of variable qi q, after some algebraic transformations and regrouping the 
terms, we result in final expression for off-diagonal contribution to radiation energy loss 
of high-energy color particle a within static approximation 

\ / off— diag. 

Here, on the right-hand side the function Ai is 

/fy \3 /(-<{») (-<W\ ^ ^ 1 

A. = (f ) C.E^ (^j„,^i:yAd..I,„CA'(.))/dq^^^ (6.7) 



X 



6{va ■ k), 
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and the function A2 has a form 



A2 = - - E ' ]np E /c^kcia;a;ImrA*(A:)){(e(k,C)-v,)(e(k,C')-v4 

(6.8) 

^ (q2+^^2)2 { *^\k)(e\k, C) T^°^°(fc, g, -k, -g)e^(k, C'))}^^ " 

In the last expression *r*°-'° is an effective four-gluon vertex entered in Paper I (Eq.(I.5.6)). 
The diagrammatic interpretation of terms in Ai and A2 that easily follows from initial 
diagrams of Fig.lUtype, is presented in Fig. El Remind that two diagrams corresponding 
to radiation of soft gluon prior to one-gluon exchange and after it, are determined by the 
first term on the right-hand side of Eq. ()6.7p . We depicted the second diagram in Fig. El 
only. Note also that functions Ai and A2 are not vanishing for plasma excitations lying 
off mass-shell due to delta-function in integrands. 

From the form of graphs in Fig.j^it is evident that they represent the contact double 
Born graphs, first considered by Zakharov J2] to ensure unitarity as was already men- 
tioned in Introduction. However in our case propagators and vertices are effective (i.e. 
we take into account hard thermal loops effects) and besides additional contribution (the 
first diagram for A2 in Fig. El) connected with existence of four-gluon HTL amplitude 
*r4g = V^g + dV^^g, appears. 

Usually the contributions to radiation energy loss containing bare four-gluon vertex 
Tig, are omitted for kinematical reasons |Sl El (the absence of momentum dependence). 
However HTL-correction 5Tig have highly nontrivial momentum dependence and therefore 
in advance, it is not evident that this term can give negligible small contribution, at least 
to region of soft transfer and soft gluon bremsstrahlung. 

Let us analyze the role of contribution Ai to the theory under consideration of brems- 
strahlung. For this purpose we compare it with basic diagonal contribution ()3.9|) (more 
exactly with the first term on the right-hand side of ()3.9j) ). Setting v^? = we rewrite 
this diagonal contribution once more, considered module squared |e'Ar*|^ 



X 



(e ■ Vc,)2 ^ (e • V, 



- 2 f-^ Re {e^(k, C) T^°'^(fc, -g, -k + qTV^Ak ' (6-9) 



+ 



e^(k, C) *r'^ik, -q, -k + qyVc^k ' 



5{va ■k + Wo.-q). 
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Figure 5: The diagrammatic interpretation of terms defining 'off-diagonal' contribution to 
radiation energy loss. The dotted lines denote thermal partons absorbing virtual bremsstrahlung 
gluons. Here, q is three-dimensional vector. 
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If one defines the last expression on mass-shell of plasma oscillations using formula ()3.10p . 
then factors 111 1 

(va ■ q)2 {va ■ ky ' (v„ ■ q) {va ■ k) 

occurring in the first two terms in square brackets of integrand will not to be singular 
due to the fact that linear Landau damping process is absent in QGP. However for off 
mass-shell excitations of medium, when frequency and momentum of plasma excitations 
approach to "Cherenkov cone" 

{va ■ fc) 0, (6.10) 

these factors become singular, that results in divergence of the integral on the right-hand 
side of Eq. fl6.9|) . From comparison of the expressions ()6.9p and ()6.7p we see that these 
singularities are exactly compensated by appropriate terms on the right-hand side of 
Eq. ()(i.7|l . Thus, the only complete sum of all contributions ('diagonal' and 'off-diagonal', 
on-shell and off-shell) to radiation energy loss of energetic parton a is a finite value. 

The meaning of the second contribution A2 (Eq. ()6.8p ) is less evident. We can only say 
that following reasoning in Section 9 of Paper II, this contribution to energy loss takes 
into account such subtle collective effect as change of chromodielectric properties of hot 
QCD medium induced by soft gluon self-interaction. This effect partially can be taken 
into account by replacement of the HTL-resummed scalar propagator *A*(fc) by effective 
one *A*(A;) taking into account nonlinear effects of soft gluon self- interaction in lowest 
order (elastic rescattering of two soft transverse gluons gg ^ gg) 

*K-^\k) - *A-i*(A;) -n«*(A;). 

Here, H*^^)*(A;) represents the correction to transverse part of gluon self-energy in HTL- 
approximation, taking into account a change of dielectric properties of medium by the 
action of the processes of nonlinear interaction (transverse) plasma excitations among 
themselves. Its explicit form is similar to (II. 9. 5) for longitudinal part of correction. 

7 Bremsstrahlung of two soft gluons 

Hereafter we are concerned with higher-order processes produced soft gluon radiation. 
Here, we consider scattering process of color particle a off hard thermal parton /3, when 
bremsstrahlung of two gluons, arises. Here, we are mainly interested in question for what 
typical amplitudes the soft gluon field (or, in other words, the density of soft gluon radia- 
tion) in medium, the process of bremsstrahlung of higher order will give a contribution to 
gluon radiation intensity of the same order as the process of bremsstrahlung discussed in 
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previous Sections. As was mentioned in Section 6, this process is defined by the effective 
current 

j;^^''[Qoa,Qofs,A^'^m = |ir;^^^(v,,v^;b,Qoa,Qo/5|fc,-fci)AW"^''n^i)c?^i- (7-1) 

We take the coefficient function in integrand in the leading order in couphng constant, 
i.e. 

KZ\ v^; b, Qoa, QoM K -h) = K^l'^i^^a, v^; b| k, -k,)QlQl^. (7.2) 

For the deriving of function K'^^^'^'^iy^, V/3; b| k, —ki) we use simple and effective procedure 
suggested in Paper I and Paper II and extended for this case. The calculations shown 
that this function has a following color structure 

KZ\'>c.. b| k-k,) = {T^T'^^Y^K,^, (v„, v^; b| -A;i) + (T«^T'^)XimK, b| -k,, k) 

The right-hand side of this expression is automatically symmetric with respect to per- 
mutation of external soft gluon legs: k ^ —ki, /i ^ /ii, a ^ Oi. The requirement of 
symmetry over permutation of hard lines {a ^ (3, c ^ b) leads to additional condition 
for partial coefficient function -R'^/ii : 

k, -ki) = K^^f,{vf3,Va; -b| - ki,k). 

The explicit form of this function, and also graphic interpretation of different terms is 
given in Appendix A. 

Let us substitute an expression for an effective current ()7.1|1 . ()7.2|1 into formula ()3.1|1 
defining total energy of radiation field in single scattering. Performing average over color 
charges, we obtain 

W2,{h) = (^^j i^^^^^\Jdkduulm{K;';^^^^^^^ (7.3) 

y<K;f;^>'\^r^,^rp■Mk-k,)){A*^'^^^^^{k^)A'~'^<^'^^ 

For the conditions of stationary and homogeneity of hot QCD plasma we have for corre- 
lation function of the soft boson excitations 

^^*(o)aiMi(^^)^(o)<K(^'j) = 5"i</^i^'i(A;i)5(A;i - k[). 
One defines the spectral density I^^{k) in the form of an expansion 

where in turn the scalar functions and are taken in the form of the quasiparticle 
approximation 

I'^' = l\.'5{u-J^') + I':i5{u + J^') (7.4) 
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In the last line we pass from functions /^'' to the number densities of soft plasma oscil- 
lations A^^''. In Eq. ()7.4p the term containing delta-function 6{uj + w^'') only will define 
contribution to energy of radiation field ()7.3p connected with bremsstrahlung of two glu- 
ons. The term containing 5(to' — c<j^'') gives contribution associated with process extending 
nonlinear Landau damping process (Section 3, Paper II) to the case of two hard color 
partons. Therefore hereafter we drop the last contribution and focus our attention on 
pure bremsstrahlung. Note that to take into account effects of weakly inhomogeneous 
and weakly non-stationary state of QGP induced by the effects of the nonlinear inter- 
action of waves and particles, it is sufficient in the first approximation to assume that 
functions A^^' slowly dependent on x = (t, x). In this case the functions obey proper 
kinetic equations, one of which will be given below. 

For simplicity we restrict our consideration to bremsstrahlung of real transverse gluons. 
With allowance for above-mentioned in ()7.3p we make following substitution 

I'^'Hh) ^ (||i^)iVik,5K + <)^E^e-(ki,Ci)e^Hki,Ci), 

C=l,2 

lmCA-'\k)) ^ -vrsignH - 4). 

Here, we have taken into account representation ()3.5p . Performing trivial integration over 
dk[, duj and dui and replacing variable ki — ki (cj^^ — > — cu^J we derive instead of 
(HI 

X I ir-^^--(v„, v^; b| k, A;i)e^(k, C) e-(ki, Ci)| L-sheii- (7-5) 

Let us define an expression for gluon radiation intensity. For this purpose we write partial 
coefficient function K^j^^-^ in the form of integral representation 

^MMi(va,V/3;b| k, ki) =Je''^'^[ICf,f,^{va,V/3\ k, ki, q)]qg=^^.^6{ujl+ujl-Vf3-q-Vai'k+ki-q))dq. 

The function /C^^^ is defined by integrand in Eq. (A.l). Subsequent reasoning are com- 
pletely similar to those in Section 3. Performing integration over impact parameter with 
()3.8j) and allowing for the following equalities 
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we obtain expressions for intensity radiation in this case 

/3=« 3 g ^ 



dA 



J2 /rfkrfki^iv*^ 



CCi=l,2 



2ui ]\ 2ujL 



X 



[j ^lf\M ^) /^/^q ^(^^ + 4. - V, ■ q - V. ■ (k + k, - q)) (7.7) 



x| I /C^^Hv., v^l fc, fci, g)e^(k, Oe^Hki, Ci)l ' + I /C^^^(va, v^l fci, fc, g)e^Hki, Ci)^^ 

+Re [( (v,, v^l fc, fci, g)e^(k, C)e^^ (ki, Ci))Yr^'(v„, v^| k, q)e'^ (ki, Ci)e^(k, ())]] 

Here, the integrand can be made completely symmetric over permutation k ^ ki if we 
replace: c^k^ki ~^ (^k^ki + '^ki^k)/2- Similar to one can define probability of 

two-gluon bremsstrahlung p^(k, ki; q) if radiation intensity is written as follows 



Comparing the last expression with ()7.7|) it is not difficult to obtain an explicit form of 
required probability. 

In limit static case of target (v/3 = 0) an expression ()7.7j] defines 'net' energy loss of 
projectile a. Setting = in Eq. (A.l) we obtain considerable simpler expression for 
partial coefficient function /C**^ entering into integrand ()7.7|) 



+ 



< *-DcuAq - h) T^'^^°(g ~h,h, -q) 



{7.1 



A' 



Va- (q- h 

+ *r''\k, fci -q,-k-ki + qyVcuuiq ~ k,) - k^, k^ -qyVcxx{k + k^ ~ q)v, 

+ *r'^\k, ki, -k-kiYVcxx{k+ki) ( ^+ T^'°'^(A;+A;i, -g, -k-k^+qyVcUk+ki-qVa 



*r-o^^ {k, -k-k, + q, -q, k^yVcUk + h- qK 



90=0 



Diagrammatic interpretation of different terms is presented in Fig.lHl As in previous cases, 
here we have not presented diagrams of the bremsstrahlung of two gluons a prior to one- 
gluon exchange. In semiclassical approximation all these process are automatically taken 
into account on the right-hand side of Eq. ()7.8|1 . 
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Figure 6: The processes of bremsstrahlung of two soft gluons in static potential model. 

The expression ()7.7|) in region of soft momentum transfer (|q| < gT) is distinct in 
value order from the expression ()3.9|) associated with bremsstrahlung of one soft gluon by 
factor 

Here, we take into account that for soft gluon o;^, |k| ~ gT. For low density of soft gluon 
radiation in medium, corresponding to the level of thermal fluctuations at soft scale gT 
(Eq. (1.6.5)), when 

\A,{X)\^^T, (7.10) 

from estimation ()7.9j) it follows that bremsstrahlung of two soft gluons is suppressed by 
more power g in comparison with bremsstrahlung of one soft gluon. 

In opposite case for the field as strong as allowed (Eq. (1.6.4)) 

\A,{X)\^T, AT*-! (7.11) 

we see from fl7.9|l that considered process (and also generally speaking all other higher 
processes of bremsstrahlung of arbitrary number of soft gluons) becomes of the same order 
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in coupling constant. Here, the problem of resummation of all relevant contributions to 
soft gluon radiation intensity arises. As mentioned above in the case of slowly inhomoge- 
neous and nonstationarity QGP, an explicit form of the number density is determined 
by solution of kinetic equation of Boltzmann type 

^ + vL ■ ^ = -NiT,[Ni] + (1 + Ni)rM], (7.12) 

where = dujl-fdh is a group velocity of the transverse oscillations. The functional 
dependence is denoted by argument of a function in square brackets. 

A generalized decay rate and regenerating rate Fi are written in the form of func- 
tional expansion in powers of the number density 



UK] = E ri"^[^k], mi] = E r!"^[^k]- (r.is) 

re=0 n=0 



Here, 



ri"^ [A^^] = ^ E / /p.- q/ «^P., P. (k, ki, . . . , k„ ■ q) AT*^ ... AT, 



ri"^[A^k]= E /7^/p./^P.,P.(k,ki,...,k.;q)(l + Ar*J...(l + Ar*J 



P=g, 1 g ' 



dq -A- dki 



^(-E'pc + Ep^ - Ep^+q-ktot - -^P,3-q - ^tot) X3 n 



(27r)3 l\ (27r)= 



where 



n n 



Stot = ^^k + E ^k,5 ktot = k + ^ ki 

1=1 i=l 

and p^(k, ki, . . . , k„; q) is a probability of bremsstrahlung of (n+ 1) soft gluons. The 
particular expressions for the cases n = and n = 1 are given by Eqs. ()H.12|1 and ()7.7j) . 
correspondingly. Note that, the expressions for generalized decay and regenerating rates 
were simplified in the sense that a contributions of more complicated 'mixed' type being 
combination of the scattering processes of Compton type and pure bremsstrahlung, are 
not taken into account. 

In conditions of ()7.1()|1 each subsequent term in the functional expansions ()7.1|1 is 
suppressed by more power g, and here we can only restrict ourselves to the first leading 
term. By virtue of the fact that T^^\ are not A^^ dependent, the Boltzmann equation ()7.12|) 
will be linear with respect to the number density. For highly excited system, Eq. ()7.1H) . 
all terms in the expansions ()7.13|1 become of the same order in magnitude, and kinetic 
equation ()7.12|1 is completely nonlinear one. 
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8 Bremsstrahlung of two soft gluons (continuation) 



Let us analyze the expression for radiation intensity ()7.7j) in more detail. For this purpose 
we introduce new partial coefficient functions 

IC^^^''' v^l k, k,, q) = ^ [k'^' (v„, v^I k, fci, q) + /C^i^(v„, v^| k^, k, g)) , (8.1) 

/C(^)**^ (v„, v^l k, ki, q) = ^ (/C"i (v,, v^l k, fci, q) - lC'\w^, v^| k^, k, q) 

the former of which is symmetric function and the second one is antisymmetric relative to 
permutation: k ^ ki, i ^ ii. In terms of the functions ()8.ip an expression within curly 
brackets in integrand of equation ()7.7|) is rewritten in the form (for brevity we suppress 
arguments of functions) 



/C^^^eVM2 + l/C'^V^eN^ + Re 



As we see from the right-hand side of this expression an interference term is vanish 
when start using new functions, and thus the probability of bremsstrahlung of two gluons 
decomposes into two completely independent parts. This points to the fact that the 
process can proceed through two physically independent channels determined by parity 
of final state of two gluon's system with respect to permutation of external soft-gluon legs. 
The functions /C'^'^)**^e*e*/ and /C'^'^^^^^e'e*/ define a matrix elements of bremsstrahlung of 
two gluons being in even state and in odd state^°, correspondingly. 

Note that going to functions ()8.1|) having more clear physical meaning can be al- 
ready performed at the level of initial effective current ()7.1|) (with approximation ()7.2|) 
for integrand) by rewriting coefficient function K'l^^^^'^ in the form 

KZtiy.^^pM k-k,) = {T'^,T'^^rK%\{y^,^,p-h\ -A;i) + [T%n''^irif)(v,,v^;b| -k,,k), 
where 

^i5i^Hva,V/3;b| fc,-A;i) = ^ (^/C^^^ (v«, V/3; b| fc,-fci) ± /^/.^^(vc,, V/3; b| - ki,k)^ 

are symmetrized and anti-symmetrized partial coefficient functions in 'impact parameter 
representation'. A purely "Abelian" part of the bremsstrahlung of two gluons is defined 
by a term {T« T«i}''^ir(2 (v,,v^;b| k-k,). 

One would expect that after separation of 'elastic' factor from /C^^j (v^, v/j| k,ki,q) 
(i.e. the factor defining a scattering process without radiation) at least symmetric function 

^"The classification of states of two plioton's system in tlieir c.m.s. is given in [JJ]. We assume tliat 
tliis classification is true for system of two gluons also. 
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+ 



^^ffji (^o' ^/sl ^' ^1' ^) "^^^ tie presented in factorized form - in the form of product of two 
functions separately depending on momenta k and ki. Such factorization takes place 
in the case of QED interaction in the semiclassical limit [ZZl- Further for simplicity we 
restrict our consideration to static case V/3 = (go = 0). Making use explicit expression 
()7.8p . we write out a matrix element /C'^'^^^^e^e^^ separating out 'elastic' factor 

/C(^)"^(v„, 0| k, fci, g)e'(k, Oe^^ki, Ci) = ,\ , ^ M^'^''^ (v,, 0| k, k,, g)e^(k, C)e*Hki, Ci), 

where 

M^'^^'^w^MkMA) = 7 Tsr—rr (8.2) 

[Va ■ k){vc, ■ ki) 

+ I < *Vcuu(q - ki) - k^, k,, -q) + {k ^ k,, z ^ 

T^-^^A;, k, -q,~k-k^ + qyVcUq - k,) T^'^^°(g - ki, k^, -qyVcxx{k + k^- q)v^ 

+ {k ^ ki, i ^ ii) 

+ ( *r'^\k, k^,~k- k,) + *r^'\k,, k, -k - k^)) *Vcxx{k + k,) 

^ { 7^ + + -q, -k-k, + qyVcUk + h- q)v'A 

- |*r^°^H^, -k-k^ + q, -q, k,) + *r^'^°^(A;i, -k - k^ + q, -q, k)yVcuu{k + k, - 

On writing this expression we use energy-momentum conservation law in the semiclassical 
limit: Va ■ {k + ki — q) = 0. The fourth term on the right-hand side of Eq. ()8.2|1 vanishes 
owing to property 

*r''^{k, ki,-k- ki) = - *r''\ki, k, -k - fci). 

The last term on the right-hand side of Eq. ()8.2j) contains contributions associated with 
effective four-gluon vertex *T4^g = T^g + ST^g. If we drop four-gluon HTL-correction ST4g, 
then this term assumes the following form in the static limit 

-2S-^j^±^^,*AKk + k,-q). 
(k + ki - q)^ 

Thus a term containing bare four-gluon vertices T^g will be proportional to longitudinal 
part of HTL-resummed gluon propagator. 

From Eq. ()8.2p we see that even if we drop the contributions containing longitudinal 
part of gluon propagators (in particular contributions contained four-gluon vertex T4^g) 
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and neglect by HTL correction 6T4g, then we still can't present symmetrized function 
A4^'^-*"^(vq,, 0| k,ki,q) in the factorized form 

A^(^)**^ (v„, 0| k, k,, q) = M'{v^, 0| k, q)W{w^, 0| fci, g), (8.3) 

where 

M'iy^, 0| fc, g) = + < *VcuM - k) T'^'*°(g - k, k, -q) 

[Va ■ k) 

is function on which the probability of bremsstrahlung of soft transverse gluon (Eqs. ()4.H) . 
()4.3p ) is determined. Let us suppose that an equality ()8.3|) is true for the time being and 
even moreover a similar factorization is true for radiation process of arbitrary number of 
soft gluons, i.e. 

A^(^)^^i---*"(v,, 0| fc, fci, . . . , g) = A<*(v„, 0| fc, g)A^^Hva, 0| fci, g) . . . A^^"(v„, 0| ^ 

(8.4) 

n = 1,2,.... 

In this case it is easy to show that an expression for total energy loss of high-energy parton 
induced by bremsstrahlung of arbitrary number of soft gluons can be presented in highly 
compact form 

i dw^!^{h) \ _^i dwgKh) \ 



X exp 



/b n=i\ dt 1^ 
^ ( Ci"'c^.\ f dcj fdr 

5 ^^T.J dke^^^-'^^ui (^1^) I >l^(v«, 0| k, g)eXk, C)| ^ 
^ CaEJ dki e^^^-'^^^K^ I (v., 0| fc„ g)e- (k^, Ci) | U 



X 



90=0 



In the second line of this equation we separate out an 'elastic' part. In the deriving of 
Eq. ()8.5|1 we making use the fact that the momentum-energy constraint can be factorized 
through the Fourier representation 



dr 



/(IT ' . 

_e^(v.-q)vK-fc)r-Qg.K.fcOr_ 

— oo 

At the soft momentum scale for low-excited QGP (an estimation ()7.1()|l ) exponential 
factor on the right-hand side of Eq. ()8.5|) is of order one and thus the energy loss in the 
leading order is defined by bremsstrahlung of one gluon. In opposite case high-excited 
QGP (an estimation ()7.11|1 ) argument of exponential function is not perturbatively small, 
and therefore in this case we can expect essential increase of radiation energy loss. 
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Unfortunately for QCD interaction even within semiclassical approximation and for 
"Abelian" part of soft gluon emission factorizations ()8.3|) . fl8.4p don't take place, and 
therefore a simple and descriptive formula ()8.5|1 for energy loss here is not correct. This 
points to the fact that multiple soft gluon radiation is not in general Poisson process and 
in this way the Poisson approximation for multiple soft gluon radiation used in a number 
of works (e.g. jTHlllIlllSl), is not correct. 

Let us derive an expression for symmetrized part of matrix element of two soft-gluon 
bremsstrahlung /C*^'^)**^(vq, 0| fc, fci, g)e*(k, (^)e*^(ki, Ci) in high-frequency and small-angle 
approximations, similar to expression A^*(vq,0| A;,g)e*(k, ^) defined in Section 4, i.e. in 
conditions when we can neglect by HTL-corrections to vertices of interaction and longi- 
tudinal part of gluon propagators on the right-hand side of Eq. ()8.2j) . 

The Coulomb factor in this approximation reads 

1 1 



ql + + — - + — 



where 



^ ^ "^UJl _ _ UJl 

" (ki + m2 + a;2M2) ' ''^^ ~ {kj^ + + xjM^)' ^ ' 

are formation times. Furthermore the functions (e ■ \oi)/{va ■ k) and (ei ■ \a)/{va ■ ki) 
entering into the first term on the right-hand side of Eq. ()8.2|1 is approximated by an 
expression of ()4.7j) type. The second term^^ is approximated by a similar way as it was 
made for term ()4.9|) . where except usual requirements in the form 

iUTf > 1, UJlTif > 1 

it is necessary to consider that inequalities 

LUiTf > 1, iUTif > 1 

are also fulfilled. Uniquely, the approximation of the transverse scalar propagator *A*(g — 
k) requires more accuracy. Instead of approximation ()4.10|) in this case we will have 

(k-q)2 ~ (k-q)l + l L+— ) ~k2 + (k-q)i + m2 + x2M2 + (8.7) 

and hence expression ()4.11|) must be replaced by 

*A*(g-fc)^ . (8.8) 

(k - q)l + m2 + a;2M2 + i^—j (k^^ + + xjM^) 



^^Hy terms hereafter we meant all expression standing within curly, square or round brackets. 
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The third term on the right-hand side of Eq. ()8.2|) is more comphcated for analysis. 
Neglecting by HTL-correction to bare three-gluon vertices and longitudinal part of gluon 
propagators, contracting with e*(k, ^)e*^(ki, (^i) and taking into account the equalities 
e ■ k = ei ■ ki = 0, we have an exact expression initial for analysis 

2ui *A*(g - h) *A\k + h- g)|-2(e ■ (ki - q))(ei • v,) + (ei • (2k - q))(e ■ v,) 



-(vc,-(k-ki+q))(e-ei)+ 



+ v. 



k + ki - q 

2 1^1 • q) 



(k2-(ki-q)2)(e-e0-(e-(ki-q))(ei-q 
e-vJ((2k+ki-q).(ki-q))-(a; + ^i)(e-(ki-q))(l-^^-p^-— ^) 



"(ki-q)2 

In the small-angle approximation next expressions hold 

e ■ k ~ e_L • (ki — — kj_) , ei • k ~ ei_L ■ (k — ^ki_L 

^ X ^ ^ Xi 

e • q ~ e_L ■ q±, ei ■ q ~ ■ q^, e ■ ei ~ ■ ei_L, 

(2k + ki - q) • (ki - q) ~ uji{2uj + c^i), 

etc. Let us consider in more detail an approximation of coefficient before scalar product 
(e ■ ei), that is equals to 

- v„ ■ (k - ki + q) + ^^l^^^^^^y (k' - (ki - q)'). (8.9) 

By the use of approximation ()8.7p it is easy to see that the expression within round 
brackets in the second term of Eq. ()8.9|) can be approximated as 



k^ - (ki - q)' ~ - (ki - q)l + m' + xfM^ +\Z±W^ + rn' + x'M 



X 



= *A-'\q-k,) + {u^-ul). 

Futhermore for denominator in the second term of Eq. ()8.9|1 we have 

1 1 vl ^(k + ki-q) 



- v„ 



(k + ki-q)2 {uj + uj,y/vl + {k + k,-q)l {uj + uiy " {uj + 
and such now we can write instead of Eq. ()8.9|) 

-v^ik-k^)^^-v^q^^+vliuj-u;,yvt + ~ + *A-'\q-h). (8.10) 

Here, in two last terms within accepted accuracy we can set Va — 1- In high-frequency 
approximation the third term in Eq. ()8.in|l is approximated as 

V^[U -UJi) c:^v^{k-ki)u + 



2oj 2uj\ 
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Let us subtitute this expression into Eq. ()8.10|) . Taking into account estimations 

M2 / 1 1 \ 



1 - fa 



Vaq\\ 



and definition of inverse propagator *A^^*(g — ki), it is easy to see tliat all terms con- 
taining m?g and in Eq. (jS.lOp exactly cancel, and such we obtain the following simple 
approximation of coefficient ()8.9|1 



2 



ki - (ki - q) 



UJ — UJi 
UJ + UJi 



(k + ki-q) 



t^ik^ - t^(ki - q)^ - {uj - t^i)(k_L ■ (ki - q)_L) 



Finally we write out an approximation for transverse scalar propagator *A*(/i; + fci — g): 



(k + ki - q)^ + m? + {x + XifW^ ' 



(8.11) 



Taking into account above mentioned we derive complete expression for matrix element 
in a small-angle approximation 



e± ■ kj 



90=0 



f8.12) 



ei_L ■ ki_L 



.)*^tf (e±- (k-q)x)(eix-kix) , 



k5_ + m2 + a;2M2 k^_L + + xfM^ 

ei± ■ (ki - q)±)(e^ ■ kj 



+ {*A\q-k,) 



+ 



Xi , 
X 



(e_L ■ (ki - -ik - q)^) (ei^ ■ (ki - q)_L) + (ei_L ■ (ki - -^k - q)^) (e_L ■ kj 
e_L • ei^)(xik^ - x(ki - q)^ - {x - Xi)(k_L ■ (ki - q)_L) 



Xi 



{x + XiY 

+ (k_L ^ ki_L, e_L ^ ei_L, x ^ Xi)\ *A\k + ki - q) 



Here, the propagators *A*(g — k), *A\q — ki) and *A*(/c + fci — g) are defined by Eqs. ()8.8|1 
and ()8.1H) correspondingly. By similar way matrix element /C*^"^)**^e*e*^ is calculated. Here, 
a new contribution defined by fourth term on the right-hand side of Eq. ()8.2p is appeared, 
where now a difference 



*r'^\k, ki,-k- ki) - *r''\ki, k, -k -ki) = 2 *r'^\k, ki, -k - ki). 
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will be stay instead of a sum of two three-gluon vertices. The explicit expression for this 
matrix element in high-energy and small-angle approximations is written out in Appendix 
A. 

Matrix elements ()8.12j) and (A. 2) can be used in particular in simulation of on-shell 
parton cascade at the build of the elliptic flow at RHIC. They enables the inelastic gg ^ 
9999 pQCD processes to be included in the consideration. Up to now only elastic gg ^ gg 
[7^ and the inelastic gg ^ ggg jHOl |HI] gluon interactions are included in simulation. 



9 Soft gluon bremsstrahlung in the case of two-scatte- 
ring thermal partons 

In this Section we extend consideration of the bremsstrahlung process to the case of 
scattering of a high-energy incident parton a off two hard thermal partons f3i and /?2 
moving with velocities V/j^ and Vjj^, correspondingly. As initial currents we choose the 
expression in the momentum space 

jSa'ik) = ^ QLv'Jiv^ ■ k), f^lZik) = ^ Ql,,y,,J{v,,,, ■ A:)e^^-°^-. 

Here, xio, X20 are coordinates of partons (3i and (32 at time t = to. In this case the energy 
of radiation field is defined (instead of (jSH}) by equation 

W{xio, X20) = -{2n)'^JdkduJdQoadQop^dQo(32 ^^Im {JQl^*''{k, xio, ^2q)*V^ {k)j''Ql^''{k, xio, X20)) 

(9.1) 

and gluon radiation intensity - by equation 

fpi3^d.Pl3i /p/sa'^Pfe dkdqidq2 
(27r)3 {2^^ (27r)9 

X 5{uj{k) - ■ qi - ■ qa - (v„ ■ (k - qi - qa)). 

Here, Wp„,p^^,p^^ is a probability of bremsstrahlung in the case of two-scattering hard 
thermal partons. The standard calculations result in following expression for an effective 
current generating this process 



X= / ^(k)wp,,p^^,p^^ (k; qi, qs) ,o,nq (9-2) 
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ja^(fc,Xio,X2o) = ^ [(T'^^T'^^)'^^Qj;„g^;3,Qfe.^M(va,v^„v^,;xio,X2o| A;) 

+ {T'-T'^r' Ql^Q%,^Q%K,{^,^, v^,; X20, xio| k) ] , (9.3) 
where partial coefficient function /('^(vq,, V/jj; xio, X20I k) is 

^/.(va, V;3^, v^2;xio,X2o| k) = I ^/^i^^/^i^ *Vc''\qi - /c)/C^'(va, 1 gi - k, -^2) 

J ■ {qi — K) 
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+ 



iv(3, ■q2)ivp, ■ (k-qi)) 



-2 



'2M 



ivi32 ■qi)iv^, ■ ik-q2)) 



+ 



Va ■ (qi + q2] 



ivp^^*Vc iqi)v(32u') {vp^x Vc (k-qi- q2)va\') 

"V(^\qi + g2)/C^'(v/3i, V/32I gi + g2,g2) 



{Va ■ {qi + q2)){Va " ^2) 

-*'^^lux{k,-ql-q2,-k + ql + q2) *^^c'''(?i + 92 ) ^i/' ( V/j^ , V/j^ | gi + ^2 , ^2 ) *'D^^\ k-qi-q2)vaX' 
- *^c'i<li + q2)v(i,u' *Vc^'{k -qi- q2)ICx'{vf3^,Vf^^ \k-q1-q2, ^2) 
+ *'i^iMuXaik, -k + qi + q2, -92, -gi) *V''(i'\k - qi - q2)vau' *'D^^'iqi)vp^x' *T^c''iQ2)vp2a'^ 

X e^^-''^°e'''^-''^°5(t;„ ■ (fc - gi - g2))5(^ft ■ gi)5(t^A ■ ^2) rfgirfg2. 

For the reason of large number of terms on the right-hand side their complete dia- 
grammatic interpretation is omitted. Here, we discuss at qualitative level a contribu- 
tion defining transition bremsstrahlung induced by four-gluon HTL amplitude (Fig.Ej): 
*T4g = r^g + Sr^g, whcre r^g is bare four-gluon vertex, and ST^g is HTL correction. In 



a 




a 



Figure 7: The transition bremsstrahlung induced by ST4g vertex HTL correction. 

the case of bare four-gluon vertex L4g this contribution is usually dropped for kinematic 
reason associated with absence of momentum dependence of L4g jSJ |3] . The momentum 
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dependence of ST^g on incoming lines is very nontrivial, and here we can expect an ap- 
preciable contribution of new mechanism of bremsstrahlung to gluon radiation intensity. 
If we restrict our consideration to this case, then under partial coefficient function in 
Eq. we mean expression 

^/.(Va, v^,, v^2;xio,X2o| k) = J /C^ ( v^, , v^^ , v^^ | k, qi, q2)5iva-{k-qi-q2))5{vf3^-qi)5{vi3^-q2) 



X e*''^° ■ e*^2° ■ ^2 dqwdq2odqi dq2 , 



(9.4) 



where 



^/.(va, v^^l k,qi,q2) = ST^^Xa{k, -gi, -q2, -k + qi + ^2) *'Dc {qi)vp^^*Vc {q2)v^, 



X *p--'(A;-gi-g2K.'. 



(9.5) 



Now we substitute an effective current ()9.3p into Eq. ()9.1|) and average over initial 
values of color charges, taking into account a relation for color algebra ()7.6|) . Keeping 
in propagator *'D'^{k) a transverse part only with allowance for ()3.5|1 and going over to 
mass-shell of radiated transverse gluon (Eq. (jH.lOj) ). after simple algebraic transformations 
we lead to the expression for radiation field energy 



iy(Xio,X2o) 



(2vr) 



'^(«)r^{/3i)^{/32)^ 
^2 ^2 ^2 

. 4 , 



t /z,(k)' 



C=l,2' 



2ul 



(9.6) 



X <j I i^*(v„, v^^; xio, X20I A;)e'(k, C)| ^ + I K\-v^, v^^, v^^; X20, xio| fc)e*(k, C)| ' 
( K'(y^, v/3j, v^^; Xio, X20I /c)e*(k, C)) (k^ (v„, V/?^; X20, Xio| fc)e^(k, Q 



+ Re 

Let us consider module squared li^'^e*^ in integrand. Substituting an expression ()9.4p 
instead of and performing integration over dqiodq2odq[Qdq2Q, we obtain 



K" ( Vc, , , ; xio , X20 1 A;)e' (k, C) r 



/C'(v«, v^,, 1 /c, gi, g2)e*(k, C) 



(9.7) 



X 



,0 _ , ^(^k-V/3i-qi-V;32-q2-Va-(k-qi-q2)) 

<?l,2-'^/3l,2'1l,2 



/C-'(v«,V/3i,v^2| A;,gJ,g2)e^(k,C) 

X 5((v, - v^J-(qi - q'l) + (v, - v^J-(q2 - q;,))e*(^^-<) ''^°e^('i^-'i^) ''^°rfqif/q2f/qWq^. 

Furthermore we analyze the expression in the last line of Eq. ()9.7j) . We present vector xio 
in the form of an expansion in terms of vectors system: the vector along relative velocity 
Vq — and the vector orthogonal to Vq, — V/j^, i.e. xio = (bi = (xio)_l, (xio)||). In this 
case a scalar product in the argument of the first exponential can be presented as 



(qi - q'l) • xio = (qi - qi)||(xio)|| + (qi - qi)± ' bi. 
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The vector bi is impact parameter of particle a relative to particle In a similar way 
we present vector X20 in the form of an expansion in terms of basis associated with vector 
of relative velocity Vq, — ^p,^: X20 = (b2 = (x2o)±, (x2o)||)- Then a scalar product in the 
argument of second exponential can be written in the form 

(q2 - Qs) • X20 = (q2 - q2)||(x2o)j| + (q2 - q2)± • W 

In this case vector b2 is impact parameter of particle a relative to particle (52- By virtue 
of the fact that vectors v^, v^^ and v^^ are fixed, such representation of scalar products 
is uniquely determined. 

With allowance for above-mentioned the following equality will be hold 

5((v. - v^J-(qi - q'l) + (v, - V;3j-(q2 - q^))e^('i-<)--e^('i-^S^)-- 



(5(1 v« - v^J(qi - q'Jii + I - v^2|(q2 -q2)|j)e 



Kqi-qi)i''^ip''(i2-q2)±-b2 



(9J 



exp 



%l 



V 



I Vq - V^J 

where in the last line we have 

1 1 



iq2 - q2j 



+ exp 



-il 



V 



iqi - qi; 



+ 



V/3iKX2oj|| - I V„ - V/32|(XioJ 



V 



After integrating Eq. ()9.8p over rfbi, dh^ and dZ, we derive 



^ 5^'H(qi - q'i)±)5((qi - q'i)ii)^^'n(q2 - q'2)±)'5((q2 - q'2)ii), 

that enables us to perform integration in Eq. ()9.7|) over c/q'^c/qg. Now we write an expres- 
sion for the gluon radiation intensity for the case of two nonrelativisic scatterers 



/dw(bi,b2,;) 



i|| ».(b„ b„ V 



bi, ba I 



(27r 



(2vr) 



Substituting ()9.(jj) in this expression^^ , taking into account ()9.7j) and ()9.8j) . we obtain a 
final expression for gluon radiation intensity 



/rfl^(bi,b2,/)' 



a. 



dt 



bi, b2, / 



IT" 



/^(a)^(/3i) x-y(/32) ■ 
/ ^2 ^2 ^2 



JpIJ\p,,\ 2^?')(/p?2/| 



d/ 



d\Pf32\ 
P02\ 27r2 



^One can slightly extend an expression for radiation intensity, if replace 



Pf<i 



(27r)3 



Pt<2 



(27r)3 



Pft'P,32 (2,^)3 (27r) 



where /p^^^^p^^ two-particle distribution function. This enables us to take into account more subtle 
effects connected with a possible correlation in motion of two thermal partons (3i and ^2- 
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X 



+Re 



4^ J 4^ II /f^ko^M jy"rfqidq25(^k-V/3rqi-V/32-q2-v„-(k-qi-q2)) 

X 1 1 /CXv«, v^,, 1 g2)e^(k, 01^ + 1 /CXv«, v^,, v^j, | A;, ga, gi)e*(k, C) | ^ (9.9) 
(/C*(va,v^,,v^J A;,gi,g2)e*(k,C)) (/C'(v«, V/j^, v^, | A;, ga, gi)e*(k, C) 



Furthermore we restrict our consideration to limit of static color centers, i.e. we 
set = = 0- (III this case variable / in Eq. ()9.8|) coincides with the longitudinal 
separation between two color static center Pi and ^2-) Taking into account Eq. ()4.2p and 
keeping in the effective propagator *'D'^' [k — qi — q2)vaa' transverse part only, we derive 
for special case ()9.5|1 

/CXv„, 0, 0| fc, gi, g2)eXk, C) Ig0=g0=o = i 2 \ / 2 2 n *^*(^' k - Qi - q2) 



^e^(k, C)5r°°^(fc, -gi, -g2, -A; + gi + g2)t;^ 
where four-gluon HTL vertex correction is defined by [7H 173] 

5P°°^(fc, -gi, -g2 - + gi + g2)Uo=,o=o (9.10) 

_ 2 [d^^ V'-v^ 1 1 

J Att V ■ k V ■ {k — qi) + ie V ■ {k — qi — q2) + ie 

Now our problem is to show that function ()9.10|1 contains singularity similar to that of 
vertex function considered above, associated with zeros of function 6, Eq. (jS.C-ij) . The 
existence of such singularity enables us at least in principle in spirit of Section 5 to simplify 
analysis of expression for energy loss^'^ ()9.9p and perform it to a great extent by analytic. 

Let us define an explicit form of integral on the right-hand side of Eq. ()9.10|) replacing 
at the moment for convenience 

k ki, k - qi ^ k2, k - qi - q2 ^ k^. (9.11) 

In this case the analog of expansion is 

, ,w ,w = Y,n\ni + Y2n\ni + Y.nlni (9.12) 

47r [v ■ ki)[v ■ k2)[v ■ k^) 

+ Zi2(<ra^ + n>*2) + ^i3(^>3 + n\ri}^) + Z23(4?^3 + ^W^)- 



Remind that for potential model of medium the expression 1)9. 9|l defines radiation energy loss of 
energetic parton a. 
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Here, rii = [k2,k3], n2 = [k3,ki] and = [ki,k2]. Chosen basis"*^^ is very convenient by 
virtue of properties 

kg ■ = 0, if s 7^ /; ki ■ rii = ka ■ = kg ■ ng. 

The exphcit form of coefficient functions Yg, Zgi, s<l, s,l = 1, 2, 3 is defined in Appendix 
B. They have following general structure 



1 ^ 



Zsi, 



(9.13) 



where functions Y^ and Zgi are given by Eqs. (B.l) and (B.2), and determinant 



6' 



k? ki ■ k2 ki ■ kg 
ki ■ k2 kl k2 ■ kg 
ki kg k2 ■ kg 



^3 



(9.14) 



represents Gramian for three-dimensional vectors ki, k2, and kg. As known vanishing of 
this determinant gives necessary and sufficient condition of complanarity of these vectors. 
Note that condition ()5.6|) also can be written in the form of Gramian, and such in this 
case condition 6 = will yields condition of coUinearity of vectors k and q'. 



We return now to variables 1)9.111) and consider equation 

k2 



k ■ q'l k ■ q'2 

qf q'l ■ q2 
ki ■ qs q'l ■ q2 



ki-q'i -'2 



q? 







(9.15) 



where q'^ = k — qi, q2 = k — qi — q2. Delta- function in the integrand ()9.9j) enables us 
to fix only g2|| = ^^k/'^a- We choose coordinate system as depicted in Fig.|2l axis OZ is 
aligned with the velocity v^^. In chosen coordinate system we have k = (|k|, 6', </)), q'j^ = 
(Iq'il, '(9, /5i), q2_|_ = (|q2_|_|, vr/2, /92), and correspondingly integration measures are 

(ik = k2(i|k| sin 9d9d(p, dq'^ = q'^^filq'^^l sm'dd'ddj32, dq^2± = |q'2±M|q'2±M/^2- 

On the conditions k^ 7^ 0, q'^^^ 7^ equation ()9.15j] can be presented as 



X 



(sin 6* sin -i? cos (v9 — Pi) + cos 6* cos ^9)" 



+ 



'1 + X^) - (x sin 6^ cos(v9 - /32) + cos ( 



+ 



(l + X^) - ix sin 13 cos(/?i -/32) + COS -df 



(9.16) 



0. 



^*If we enter vector N = [ki,ni] + [k2,n2] + [k3,n3] orthogonal to vectors ki, k2, and ka, then this 
expansion basis should be added to tensor N''N^ . However by Jacobi identity we have N = 0. We also 
can enter identity tensor XS^Hn expansion l|9.12(l . As shown in Appendix B a coefficient of expansion X 
without loss of generality can be set equal to zero. In this sense chosen expansion basis is complete. 
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Here, x = '^a|Q2±l/'^- This equation represents appropriate extension of Eq. ()5.7|) . Note 
that angle j32 will be enter everywhere as difference with azimuth angles ip and and 
therefore it can be set equal to zero and integral is / d(32 = 27i. 

We can obtain simplest nontrivial solution of this equation if we set expressions in two 
last square brackets on the left-hand side of Eq. ()9.16|) equal to zero. The solution have a 
form 



cos(p± = — cos^T a/1 + xM, cos/?i± = — cos^9ta/1 + xM. (9.17) 

It is not difficult to see that in deciding of ()9.17p the expression in first square brackets 
vanishes. This solution is correspondent to the case already considered in Section 5, 
and therefore performed analysis of transition bremsstrahlung can be extended to this 
problem. However solution ()9.17j) represents very hard condition for vanishing of function 
6'. It corresponds to condition of coUinearity of three vectors k, q'^ and q^2- general case 
it is necessary to solve algebraic equation ()9.16p . that in notation x = cosy?, y = cos/5i, 
can be written in more clear form 



where 



[xy - ^(l-x2)(l-7/2)) + 6] + {cex + def + {c^y + d^f -3 = 0, 



sin 6 sin i}, b = cos 6 cos i), ce = tanh ^ sin 6, = tanh ^ sin 
cos 9 , cos 



, cosy , cosv ^ , / r, o\ ■ , ^ 

de = 777' "1? = TT' ^ = Mx + y^ + X^), X = smh^. 

cosh 4 cosh 4 ^ * ^ 

At the end of this Section we briefiy analyze a case of transition bremsstrahlung in 
scattering of energetic parton a off three static color centers, the process proceeding 
through five soft gluon HTL induced vertex, as depicted on Fig.|Hl In the case of radia- 
tion of transverse soft gluon, calculation of matrix element of this process is reduced to 
calculation of integral in the form 



An {v ■ ki){v ■ k2){v ■ k^){v ■ k^y 

where ki = k, k2 = k — qi, k^ = k — qi — q2 and k^^ = k — qi — q2 — q-i. According reasoning 
of Frenkel and Taylor [73] , this integral represents pure spatial part of Lorentz-covariant 
tensor of 2-rank 

J An {v ■ ki){v ■ k2){v ■ ks){v ■ k4,) 
This integral can be presented in the form of expansion in basis 

[/^^ = ^AX< + E^^K« + «), s,/ = 1,2,3,4, (9.18) 

S = l S<1 
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vp^=0 Vp^=o Vp^=o 



Figure 8: Transition bremsstrahlung induced by ST^g vertex HTL correction. 



where we enter 4- vectors 



— b A.4 , 




c r\j^ th^ th^ , 



,3 — fc l^il^i 1^2 ) 




/X!/A(T 7 1/ 7 A 7 cr 



and 5^*'^-'*°" is absolute antisymmetric tensor of 4-rank. In Appendix C calculation of 
coefficients and Bgi, s < I are given. They have following general structure analogous 



where the functions Ag and Bgi are given by expressions (C.l) and (C.2), correspondingly, 
and 



represents Gramian for 4- vectors ki, ^2, k^ and k^. Zeros of determinant ()9.20|) define sin- 
gular surface in a space of kinematic variables, on which it should be defined all functions 
in matrix element of this transition bremsstrahlung. 

If we consider medium induced bremsstrahlung of higher-order than in Fig.|Hl i.e. one 
in scattering off four, five etc. static color centers, then we can expect an absence of general 
singular factor as it occurs in coefficient functions ()5.4|1 . ()9.13|1 and ()9.19|) by virtue of 
the fact that it is impossible to define nontrivial Gramian for five and more 4-vectors. 
Therefore we can suppose that these higher processes of transition bremsstrahlung will 



to ^rn^ 




(9.19) 




(9.20) 
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be strongly suppressed in comparison with lower-order processes (i.e. containing gluon 
vertices with smaller number of soft external legs) and thus they are not of our interest. 

10 Summary and outlook 

In this work in the framework of semiclassical approximation we have proposed an ap- 
proach to calculation of radiation intensity of high-energy parton travelling through hot 
QCD medium. The deriving of some effective color currents generating bremsstrahlung 
processes of arbitrary number of soft gluons for scattering energetic parton off arbitrary 
number of thermal particles of medium is a great moment for this consideration. 

A suggested scheme for calculation of these effective currents tacitly implies by virtue 
of initial dynamical equations that hot QCD medium represents a weakly coupled sys- 
tem (^f -C 1), allowing to use perturbative approach. However theoretical analysis of 
experimental data, especially concerning with collective flow, unambiguously shows that 
properties of 'new matter' produced at RHIC experiments for moderate temperatures 
(1 ^3)Tc substantially different from the expectation of a weakly coupled colored plasma 
(see the last reviews on this subject [H ^ I8H | IM | l85j). A major revision of proposed scheme 
to calculation of radiative energy loss can be one of the consequences of such conclusion. 
For example within standard Gyulassy-Wang model multiple scattering in medium is con- 
sidered as totality of subsequent acts of pair interaction of particles, and in particular an 
influence of simultaneous collision of several (more then two) particles on formation of 
bremsstrahlung spectrum is neglected. It is clear that effects of such type can play an 
important role in the processes of gluon multiplication in rather dense scattering medium 
(e.g., liquid type), when plasma parameter F becomes more or of the order of one. The 
attempt of taking into account of such effects in dense matter was made in work [Hn] for 
the case of bremsstrahlung of fast electric charged particles. Here, diagrammatic formal- 
ism for two-particle Green functions was used. It was shown that in the case of radiation 
in sufficiently dense medium in a distant long wavelength region of spectrum coherent 
multiple scattering of particles results in (supplementary to LPM effect) suppression of 
the yield of bremsstrahlung photons. We can expect that similar phenomenon has to 
take place and for the case of dense^^ strongly coupled quark-gluon plasma (sQGP) for 
description of that a number of approaches is already proposed |HHl lEHl EOl EI] • 

^^There is some ambiguity in definition of sQGP, obeying the classical statistics, in terms of plasma 
parameter F. The ambiguity of such a kind was first noted in work ^87^ . The initial definition of F is 

r-^, (10.1) 

where n = J^b "^P ^ mean density of system and rn is Debye length, that in the case of the classical 



54 



The approach to calculation of effective currents suggested in this work (and in Papers 
I and II also) can be extended to the case of strongly coupled QGP. The use of coUisionless 
Blaizot-Iancu equation (1.3.3) correctly describing behavior of QGP in the weakly coupling 
regime and at soft momentum scale of plasma density fluctuations is one of the most 
key moments for the derivation of effective color currents. Therefore as a first step of 
construction of the radiation energy loss theory in sQGP this kinetic equation is necessary 
to be extended in a way. As one of the leading ideas representations developed in attempt 
to describe thermodynamical properties, the transport coefficients and the dynamical 
correlations of strongly coupled usual plasma, can be used. 

The extensive research of properties of usual plasma, when plasma parameter F 
changes within huge region: from slowly coupling regime F <^ 1 and up to F ~ 1000 were 
undertaken more than 20 years ago. In particular, for a strongly coupled one-component 
plasma Vieillefosse and Hansen [02] and Wallenborn and Bans have advanced theoret- 
ical estimations of the shear r] and the bulk viscosities on the basis of hydrodynamic and 
kinetic-theoretical arguments, correspondingly. The two calculations are in agreement 
shown that the transverse viscosity rj (exactly its dimensionless combination) exhibits a 
minimum as a function of F (for values F ~ 10 ^ 20) and ( is negligible compared to 
T]. The value r/ in a point of minimality is more little then a value for the weakly cou- 
pled plasma, when F <^ 1. This very important and striking observation was confirmed 
through molecular-dynamics computer simulations [HI] . The position of the minimum is, 
however, rather sensitive to the various approximations. 

Furthermore Tanaka and Ichimaru |95t l96] have developed a new kinetic theory of dy- 
namic correlations for a strongly coupled, classical one-component plasma within the gen- 
eralized viscoelastic formalism [DTJ . The strong Coulomb-coupling effects beyond random 

non-Abelian plasma reads 

rl- , (10.2) 

0=q, g, g 

fcs is a Boltzmann constant. In the theory of usual plasma strong coupling regime is defined by condition 
r > f . In the case of QED interaction the coupling is weak, and thus the condition F > 1 can be realized 
only for large value of plasma density. For QCD interaction a situation is more complicated. As it is seen 
from Eqs. H10.1|l . (|10.2(l the last condition can fulfilled in general and for a rarefied plasma, but with the 
strong coupling between thermal particles and for very dense plasma but with weak coupling. For the 
quantum non-Abelian plasma, when 

n 1 / he \^ ( knT 



\ fc^T / ' \ he 

such ambiguity does not exist. Here, the plasma parameter is more or equals one if and only if g > 1 j87) . 

To avoid the ambiguity in definition of sQGP for classical statistics we will assume the medium is dense 
and the coupling is strong. 
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phase approximation are described in the present theory through the dynamic local-field 
correction G(k, uj) introduced via linear-response relation between the external potential 
V9ext(k, 1^) and the induced density fiuctuation 5n(k, cu). Making use a fully convergent 
kinetic equation originally proposed by them, where higher-order scattering processes 
are included through the factor 1 — G(\<.,uj) in integrand of collision term, Tanaka, and 
Ichimaru have established a self-consistent scheme to calculate t] and G(k, u) simultane- 
ously. Their predictions showed a good agreement with early obtained results of different 
approaches jnaEHlini. 

Performed researches for strongly coupled usual plasma suggest that in a case of classi- 
cal sQGP we can also expect an existence of minimal value of the shear viscosity for some 
finite value of QCD plasma parameter p0.1|) and this value can be appreciably smaller 
than value following from perturbative estimations 1^ for weakly coupled QGP. It 
is represented extremely important by virtue of very small value of t] for a new QCD 
matter discovered at RHIC experiments |l()()j . The requirement of minimality of shear 
viscosity has to lead to definite fixing of generalized collision term for kinetic equation, 
that in weakly coupled limit must turn to kinetic equation (1.3.3). Having in hand such 
an equation we can further develop an extended scheme for calculations of color effective 
currents generating bremsstrahlung processes (and radiative energy losses produced by 
them) of high-energy partons in the sQGP. The exact consideration of this approach will 
be the subject of a further research. 
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Appendix A 



Here, we give complete expression for partial coefficient function ^^'^^^(vq,, V/3; b| A;, fci) 
defining bremsstrahlung of two soft gluons in scattering of energetic color particle a off 
hard thermal parton /3 

K^,^^{\a, V/3; b| k, ki) = \ — N - h)K^^^^ (v/jI g - fci, /C2) 

J lVj3 ■ [q — Ki) 

+ *T^,x{k, ki -q,-k-k^ + q) *VlY'{q - ki)K,,^,{^p\ q - k,, k{) *V^c^'{k + k^ - g)t;„A 
+ *T^^,x{k, ki, -k - ki) *V^^\k + ki)Kx'{ya, v^| /c + fci, g) 
- T^,ami(^, -k-k^ + q, -q, k,) ^V'S'lk + k, - q)vp,, *V^^\q)vpx^ 

X e"*^- 5{v^ -{k + ki- q))6{vp ■ q)dq. 

Here, the functions K^^^ (v| /c, ki) and /C^(vq, v^| k, q) are defined by Eq. ()2.10|) and Eq. ()3.7p . 
correspondingly. The graphic interpretation of various terms in this expression is pre- 
sented on Fig.ini Remind that graphs, where radiation lines from hard lines are located a 
prior one-gluon exchange, are not given here. 
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Figure 9: Two soft gluon bremsstrahlung. 
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Below we adduce a complete expression for matrix element /C^^^^^^^e^e*^ in the high- 
frequency and small-angle approximations. This matrix element defines bremsstrahlung 
of two gluons forming odd state system (in the c.m.s. of gluons) with respect to space 
inversion 

/C(^)^^^ (v„, 0| k, h, q)e\k, Oe^Hki, Ci)l .0=0 



2 M 

q± + /^D + — — + — 



e_L ■ kj 



ei_L ■ k 



+ *A*(g - k) 



(ex ■ (k-q)^)(eix -kix) , , , (ei± • (ki - q)x)(ex ■ k^) 



+ {*A\q-k,) 



kj^ + m2 + xf M2 



Xi 



*A\q - k. 



X 



e_L ■ ki k-q (ei^ ■ (ki -q)_L) + ei_L ■ ki k-q (e_L ■ k 



k^ + m2 + x2M2 

Xi 



X 



+ Xi 



e_L ■ ei_L, 



xik^ - x(ki - q)^ - (x - Xi)(k_L ■ (ki - q) 

[X + XiY L 

+ (k^ ^ ki_L, ^ ei_L, x ^ */\\k + ki - q) 



(A.2) 



*A\k + k,)(^^ + '^^ 



-(l + ^](e,..(k-£k4)(e..k. 



X 
Xi 



l + ^](e,.(k,-^k)J(e,,.k,, 



2^, 

X 



(e± • eix) r 



X + Xi 



xikl_ - a;ki_L - (x - 2;i)(k_L ■ ki^) 



2 *A*(fc + k^) *A\k + k,- g)| (e^ ' (^i - ^ ■ ((l + ^)ki - q 



--(k-^kOJ(e..((l + ^)k-q)J 



+ 



(e_L ■ ei^J 



xik^ - xk^_L - (x - xi)(k_L • kix) - ((xik - a;ki)_L ■ q^) 



where 

*A\k + k,] 



(k + ki)i + m2 



1 + 

X / 



X 



(ki + m )+ 1 + - (kf^ + m 



The functions Tf, Tif are defined by Eq. ()8.8|) and the propagators *A*(g — A;), *A*(g — ki] 
and *A*(/i; + fci — g) are defined by Eqs. ()8.8|) and (jHHH), correspondingly. 
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Appendix B 



The coefficients Yg, s = 1, 2, 3 in the expansion of ()9.12|1 are defined by contracting with 
klki- Such, for example a contracting with k\kl with use of identity 

(v ■ ki)^ = Ldf — 2uji{v ■ ki) + {v ■ ki){uJi — v • ki) 

and formulas (3.24), (3.30) and (4.15) from [Hj yields, 

= Yi= ojl[auJi + huj2 + ecus] - UiM{k2, ks) (B.l) 

+ ^ {(t^2(ni ■ n2) + t^3(ni ■ n3))M(A;2, ks) - (rii ■ n2)L(A;3) - (iii ■ n3)L(A;2)|, 
where 6' is defined by Eq. fl9.14p and 

a = a{ki,k2,ks) = ^, h = h{ki,k2,kz) = c = c(A;i, /c2, /cs) = (B.2) 



Ai 



M(/c2, k'i) ki ■ k2 ki ■ ks 
M{ki,k3) kl k2-k3 
Miki,k2) k2-k3 



ki 



A, 



A' 

kl M{k2,k3) ki-ks 
ki-k2 M{ki,k'i) k2-ks 
k^-ks M{ki,k2) kl 



A. 



kl ki-k2 M{k2,k3 

ki-k2 kj M{ki,k3 
ki-ks k2-k3 M{ki,k2 



A' 



kl kl ■ ^2 kl ■ k^ 
kl ■ k2 kl k2 ■ k^ 
kl ■ ks k2- ks 



kl 



and functions M{k,ki), L{k) are defined by Eq. (jS.Sj) . The coefficients Y2 and Y3 are 
obtained in a similar way. They can be obtained from (B.l) by cyclic permutation of 
1,2,3. 

The coefficients Zgi, s < I are defined by contracting ()9.12j) with klkj. For example, a 
contracting with fc^fcg yields an expression for Z12 

6'Zi2 = Z12 = uJiUJ2[ai^i + buJ2 + cuj-:^] - uJiM{ki, k^) - uj2M{k2, k^) + -^(^3). (B.3) 

The remaining coefficients Z13 and Z23 follow from (B.3) by cyclic permutation of 1,2,3. 

Now we add a term X5^^ to the right-hand side of ()9.12j) . In this case instead of 
relations (B.l) and (B.3) we will have 

+ k^X = 5'Zi, + (k, ■ kz)X = Zu, s<l; s, I = 1,2, 3. (B.4) 

The equation for X is derived by taking of trace ()9.12|1 



X = 



Att {v ■ ki){v ■ k2){v ■ k^) 
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aui + bu!2 + CUJ3 



(B.5) 



S=l S<1 

If we substitute (B.4) into the rightmost expression in (B.5) and collect similar terms, 
then we will have an equation for X 



3i'-(i:kM + 2E(k.-k,)(n.-n,: 

S=l S<1 



1 

6' 



S = l S<1 



. (B.6) 



Let us consider an expression in square brackets in the first term on the right-hand side 
of Eq. (B.6). By virtue of the fact that scalar production represents algebraical adjunct 
to element (ij) in Gramian ()9.14|1 . then 

k^n^ + 2^(ks ■ k/)(ns • n/) = 5' (no summation over s\) (B.7) 

for all s = 1,2,3. Therefore a factor after X exactly vanishes. Let us show that the 
following identity is true 



X 



6' 



S<1 



si 



s=l 



(B.8) 



If one substitute instead of Yg and Z^i their integral representation, then the right-hand 
side of Eq. fl6.8p can be written as 



1 rdflv 
J' ~ 



3 \ 2 

E n,(v ■ k^ 

s=l 



(B.9) 



47r {v ■ ki){v ■ k2){v ■ ks) 
Let us expand unit vector v in basis {ng}: 

3 

V = XI '^^^'^^ 
s=l 

where Ks are coefficients of expansion. Then by the fact of independence of scalar product 
ks ■ lis from s and condition = 1, we will have 

^ s=l ' ^s=l ' ^ s=\ ' 

Comparing (B.9) with (B.5) we see that equality (B.8) is really true. Incidentally its 
correctness can be shown by direct substitution of the functions (B.l), (B.3) and (B.5), 
and by use of equality^^ 

A' = (nicji + W2^2 + natus)^ - 5'. 

Such we see that without loss of generality one can set coefficient of expansion X equal 
to zero, although in particular calculations it may be more convenient to use its explicit 
value. 



This equality is generalization of A — (5— (wik2— Ci'2ki)^ defined for functions A — fc^A:| — (/ci-fc2)^, 8 
k^ki - (ki • k2)2 (Eq. dSSIl). 
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Appendix C 



The expression for coefficients Ag in the expansion ()9.18|1 is obtained by contracting with 
fc^fc^, where now s = 1, 2, 3, 4 and using properties 

ks-ni = 0, for s^l, {k^ ■ n^ f = A", 

where A" is defined by Eq. ()9.20|) . Thus for Ai we derive equation 
'dQ^ (t> ■ ki) 



A"Ai = Ai 



ai{ki ■ k2) + hi{ki ■ k^) + ci(/ci ■ k^). [C.l] 



An {v ■ k2){v ■ k^){v ■ ki) 

Here, ai = 01(^2, ^3, ^4)5 W = bi{k2, k^, k^j, Ci = Ci{k2, fcs, k^) and functions ai, bi and Ci 
are defined by Eq. (B.2). The other coefficients A2, A^ and ^4 are calculated in similar 
way. 

The coefficients Bsi are derived by contracting ()9.18|1 with k'^ki- Here, we will have, 
for example, for B12 



A"B 



12 



B 



12 



An {v ■ k-i){v ■ ki) 



M{k3,ki). 



(C.2) 



Now we add term Cg'^'^ to the right-hand side of ()9.18p . In this case instead of (C.l) 
and (C.2), we will have equations 

A"^ + CA;2 = As, A"B,i + C{K ■ h) = 3^2, s < I; s,l = 1,2, 3, 4. (C.3) 

The equation for C is defined by taking the trace ()9.18|1 and allowing for t>^ = 0. Substi- 
tuting the coefficients Ag and Bsi in obtained equation and collecting similar terms, we 
derive desired equation for coefficient C 

4 T 1 r 4 



1 

'a^ 



Y,nlAs + 2j2ins-ni)Bsi 

S=l S<1 



1 

A" 



AA"-{J2kys + ^J2iks-ki)ins-ni: 

S=l S<1 



C. (C.4) 



Let us show that expression in square brackets on the left-hand side of Eq. (C.4) 
vanishes. Substituting instead of Ag and B^i their integral representations of (C.l) and 
(C.2) type, we obtain 



Y,nlA, + 2Y,{n,-ni)B,i= / 

s=l s<l •' 



4 \ 2 

'Eusiv ■ ks) 

s=l 



dflv 

An {v ■ ki){v ■ k2){v ■ k-i){v ■ ki) 



(C.5) 



Furthermore we use the same trick as supposed in Appendix B. We present 4- vector v in 
the form of expansion over basis {n^}: 

4 

s=l 



62 



where 7]s are coefficients of expansion. Then we can read 

by virtue of condition v"^ — Q and independence of ks ■ Ug on s. Such we show that 
expression on the left-hand side of Eq. (C.4) vanishes. It can be also shown by direct 
computing with use of explicit expression for Ag and Bgi of (C.l) and (C.2) type. 

The expression in square brackets on the right-hand side of Eq. (C.4) vanishes by 
equality 

kikl -\- 2 ^(/Cs • ki){ns ■ rii) — A" (no summation over s!), 

being analog of identity (B.7). Such we show that coefficient of expansion C by virtue of 
its arbitrariness can be assumed equal to zero. 
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